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ON THE RIESZ BASIS PROPERTY OF ROOT VECTORS SYSTEM FOR 

2x2 DIRAC TYPE OPERATORS 

ANTON A. LUNYOV AND MARK M. MALAMUD 


Abstract. The paper is concerned with the Riesz basis property of a boundary value problem 
associated in [0,1] 0 with the following 2x2 Dirac type equation 

Ly =+ Q(x)y = \y, ^ ^ 2 }’ 2/=)’ (^.l) 

with a summable potential matrix Q € T^[0,1] ® and &i < 0 < 62- If 62 = —bi = 1 this 
equation is equivalent to one dimensional Dirac equation. It is proved that the system of root 
functions of a linear boundary value problem constitutes a Riesz basis in L^[0, 1] 0 provided 
that the boundary conditions are strictly regular. 

By analogy with the case of ordinary differential equations, boundary conditions are called 
strictly regular if the eigenvalues of the corresponding unperturbed {Q — 0) operator are asymp¬ 
totically simple and separated. As distinguished from the Dirac case there is no simple algebraic 
criterion of the strict regularity whenever 61 -I- 62 ^ 0. However under certain restrictions on 
coefficients of the boundary linear forms we present certain algebraic criteria of the strict regu¬ 
larity in the latter case. In particular, it is shown that regular separated boundary conditions 
are always strictly regular while periodic (antiperiodlc) boundary conditions are strictly regular 
if and only if 61 -|- 62 7^ 0. 

The proof of the main result is based on existence of triangular transformation operators for 
system (0.1). Their existence is also established here in the case of a summable Q. In the case 
of regular (but not strictly regular) boundary conditions we prove the Riesz basis property with 
parentheses. The main results are applied to establish the Riesz basis property of the dynamic 
generator of spatially non-homogenous damped Timoshenko beam model. 


Contents 


1. Introduction 2 

2. Triangular transformation operators 5 

2.1. The Banach spaces Xi and X^o 5 

2.2. Transformation operators 8 

3. Asymptotic behavior of solutions 16 

4. Regular boundary conditions 19 

5. Strictly regular boundary conditions 23 

6. The Riesz basis property of root vectors system 28 

6.1. Some auxiliary results 28 

6.2. Proof of the main result 29 

6.3. The case of general potential matrix 33 

7. Application to the Timoshenko beam model 35 


2010 Mathematics Subject Classification. 47E05, 34L40, 34L10, 35L35. 

Key words and phrases. Systems of ordinary differential equations; regular boundary conditions; transforma¬ 
tion operators; Riesz basis property; Timoshenko beam model. 



2 


ANTON A. LUNYOV AND MARK M. MALAMUD 


References 


38 


1. Introduction 

Spectral theory of non-selfadjoint boundary value problems (BVP) on a finite interval X = 
(a, b) for nth order ordinary differential equations (ODE) 

y(")+W"“'^ + - + 9n-i2/ = A"y, x€(a,5), (1.1) 

with coefficients qj G L^[a,b] takes its origin in the classical papers by Birkhoff [4, 5] and 
Tamar kin [50, 51, 52], They introduced the concept of regular boundary conditions for ODE and 
investigated the asymptotic behavior of eigenvalues and eigenfunctions of related BVP. Moreover, 
they proved that the system of root functions, i.e. eigenfunctions and associated functions, of 
the regular BVP is complete. Their results are also treated in the classical monographs (see [41, 
Section 2] and [14, Chapter 19]). 

More subtle is the question of whether the system of root functions is a Riesz basis in L^[a, 
V.P. Mikhailov [38] and G.M. Keselman [21] independently proved that the system of root func¬ 
tions of a boundary value problem for equation (1.1) forms a Riesz basis provided that the 
boundary conditions are strictly regular. Similar results are also obtained in [14, Chapter 19.4]. 
Moreover, for boundary conditions which are only regular but not strictly regular, A.A. Shka- 
likov [45, 46] proved that in the case qj G L^{a,b), j G {l,...,n — 1}, the system of root 
functions forms a Riesz basis with parentheses. Recently A.M. Minkin [39] proved that the 
converse statement is almost true. Namely, he proved that if multiplicities of the eigenvalues are 
uniformly bounded, the Riesz basis property for the system of root functions of BVP implies 
the regularity (not necessarily strict regularity) of the boundary conditions. 

Numerous papers are devoted to the completeness and Riesz basis property for the Sturm- 
Liouville operator (see the recent review [31] by A.S. Makin and the references cited therein). We 
especially mention the recent achievements for periodic (anti-periodic) Sturm-Liouville operator 
— ^ [0)^]- Namely, E. Gesztesy and V.A. Tkachenko [15, 16] for q G L^[0,vr] and 

P. Djakov and B.S. Mityagin [12] for q G 1T“^’^[0, vr] established by different methods a criterion 
for the system of root functions to contain a Riesz basis. 

In this paper we consider a special case of the following first order system of ODE 

Ly := L{Q)y:=-iB~^y'+ Q{x)y = Xy, y = col(yi,..., y„), (1.2) 

where B is a nonsingular diagonal n x n matrix with complex entries, 

B = diag(6i,62,...,M ( 1 - 3 ) 

and Q{-) =: {qjk{-))'j k=i G T^([0,1]; C”'^”) is a potential matrix. 

To obtain a BVP, we adjoin to equation (1.2) the following boundary conditions (BC) 

Cy(0) + Dy(l) = 0, C = (c,fc), D = (d.^) G (1.4) 

Moreover, in what follows we always impose the maximality condition rank(C D) = re. 

Note that, systems (1.2) form a more general object than ordinary differential equations. 
Namely, the reth-order differential equation (1.1) can be reduced to the system (1.2) with r = re 
and bj = exp {27rij/n) (see [33]). The systems (1.2) are of significant interest in some theoretical 
and practical problems. Por instance, if re = 2m, B = diag(/m ,5 —Im) and Qu = Q 22 = 0, 
the system (1.2) is equivalent to the Dirac system [26, Section VII.1], [35]. Note also that 
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equation (1.2) is used to integrate the problem of N waves arising in the nonlinear optics [42, 
Section III.4], 

With the system (1.2) one associates, in a natural way, the maximal operator Lmax = .hmax(Q) 
acting in L^([0,1];C") on the domain 

dom(L^,^) = {ye WiH[0,l];C-) : G L2([0, 1]; C*^)}. (1.5) 

We denote by Lc,d '■= Lc,d{Q) the operator associated in i^([0,1];C”) with the BVP (1.2)- 
(1.4). It is defined as the restriction of L = L{Q) to the domain 

dom(Lc,D) = {y G dom(Lmax) : Cy{0) + Dy{l) = 0}. (1.6) 

Apparently, the spectral problem (1.2)-(1.4) has first been investigated by G. D. Birkhoff 
and R. E. hanger [6]. Namely, they have extended certain previous results of Birkhoff and 
Tamarkin on non-selfadjoint BVP for ODE to the case of BVP (1.2)-(1.4). More precisely, they 
introduced the concepts of regular and strictly regular boundary conditions (1.4) and investigated 
the asymptotic behavior of eigenvalues and eigenfunctions of the corresponding operator Lc,d- 
Moreover, they proved a pointwise convergence result on spectral decompositions of the operator 
Lc,d corresponding to the BVP (1.2)-(1.4) with regular boundary conditions. 

The completeness problem of the root vectors system of general BVP (1.2)-(1.4) has hrst 
been investigated in the recent paper [34] by one of the authors and L.L. Oridoroga. In this 
paper the concept of weakly regular boundary conditions for the system (1.2) was introduced 
and the completeness of the root vectors for this class of BVP was proved. Por the Dirac type 
system [B = B*) the concept of weakly regular boundary conditions (1.4) coincides with that 
of regular ones and reads as follows: 

det(CP+ + DP.) / 0 and det(C'P_ + DP+) ^ 0. (1.7) 

Here P+ and P. are the spectral projections onto ’’positive” and ’’negative” parts of the spectrum 
of B = B*, respectively. In the recent papers [30, 29] the completeness of root vectors was 
established for certain classes of non-regular and even degenerated boundary conditions under 
certain algebraic assumptions on the boundary values Q(0), Q{1), of the matrix Q{-). 

Eurther, if Dirac type operator Lq^d is dissipative, then regularity of conditions (1.4) is 
equivalent to the first of conditions (1.7) only. It is proved in [27] that the resolvent (Tc,d — A)“^ 
of any complete dissipative Dirac type operator Lc,Z) admits the spectral synthesis. In particular, 
the latter happens if the first of conditions (1.7) holds. 

Einally, in [29, 30] it was established the Riesz basis property with parentheses for system (1.2) 
subject to various classes of boundary conditions with a potential Q G L°°([0,1];C”^”). In [40] 
a stronger result was obtained for the Dirichlet BVP for 2m x 2m Dirac equation (n = 2m, 
B = diag(Im, —dm)) with a potential matrix Q G T^([0,1]; 

In this paper we investigate the Riesz basis property for 2x2 Dirac type system 

-iB~^y'+ Q{x)y = \y, y = col(yi, 1 / 2 ), a:G[0,1], (1.8) 

subject to regular and strictly regular boundary conditions (1.4). Here 

B = dmg{biM), 6 i< 0 < 62 , and Q = 1 ^ (1.9) 

Eirst we note that in this case boundary conditions (1.4) are regular, i.e. conditions (1.7) are 
valid, if and only if they are equivalent to the following conditions 

Ui{y) = yi(0) + 6^2(0) + ayi(l) = 0, U2{y) = dy2(0) + cyi(l) ^2(1) = 0, 


( 1 . 10 ) 
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with certain a,b,c,d G C satisfying ad — be ^ 0. Clearly separated, periodic, and antiperiodic 
boundary conditions are regular. 

Next we recall that regular BC (1-10) are called strictly regular, if the sequence Aq = 
of the eigenvalues of the unperturbed {Q = 0) BVP (1.8)-(1.10), is asymptotically 
separated, i.e. there exist <5 > 0 and no G N such that 

for j^k and \j\,\k\>no. 

In particular, the eigenvalues are geometrically and algebraically simple. 

For Dirac operator {B = diag(—1,1)) the strict regularity of BC reads as follows: (a — / 

—46c. 

Going over to BVP ( 1 . 8 )-( 1 . 10 ) we note that a special case of 2 x 2 Dirac operators Lc,d{Q), 
have been investigated much deeper. For instance, P. Djakov and B. Mityagin [10] imposing 
certain smoothness condition on Q proved equiconvergence of the spectral decompositions for 
2x2 Dirac equations subject to general regular boundary eonditions. 

Moreover, the Riesz basis property for 2 x 2 Dirac operators Lc,d{Q) has been investigated 
in numerous papers (see [54, 55, 18, 8 , 3, 9, 11, 12, 13] and references therein). The most 
complete result was obtained by P. Djakov and B. Mityagin in [9]. Namely, assuming that 
Q G T^[0,1] it is proved in [9] that the system of root vectors of the Dirac operator L{Q) 

with regular boundary conditions constitutes a Riesz basis with parentheses in L^[0,1] ® and 
ordinary Riesz basis provided that BC are strietly regular. Note, that non-degenerate separated 
boundary conditions are always strictly regular, hence the root vectors of the corresponding BVP 
constitute a Riesz basis [9] (see Remark 6.10 in this connection). 

The following theorem is the main result of the paper. 

Theorem 1.1. Let Lc,d{Q) be the operator associated in L^([0,1 ];C^) with the BVP (1.8)- 
(1.10) and let Qi 2 ,Q 2 i G -^^[0,1]. Assume that boundary conditions (1.10) are strictly regular. 
Then root veetors system of the operator Lc,d{Q) forms a Riesz basis in L?‘[Q, 1] (g) C^. 

While definition of the strict regularity in the case of bi —62 is rather implicit, for certain 
classes of boundary conditions it can be expressed in purely algebraic terms. For instance, if 
be = 0 and ad 0, then BC (1.10) are strictly regular whenever 6 iln|d| + 62 ln|a| 7 ^ 0. In 
particular, periodic {a = d = —1) and antiperiodic {a = d = 1) BC are strictly regular if and 
only z/ 61 + 62 7 ^ 0. Therefore Theorem 1.1 implies the following surprising result. 

Corollary 1.2. Let Q12, Q21 € A^[0,1] and 61 + 62 7^ 0. Then the system of root vectors of the 
periodic (antiperiodic) operator Lc,d{Q) forms a Riesz basis in T^[0,1] (g) C^. 

This result demonstrates substantial difference between Dirac and Dirac type operators. Note 
in this connection that a criterion for the system of root functions of the periodic (necessarily 
non-strictly regular) BVP for 2x2 Dirac equation to contain a Riesz basis (without parentheses) 
was obtained by P. Djakov and B. Mityagin in [12]. 

We also prove that the root vectors system of the operator L{Q) forms a Riesz basis with 
parentheses provided that BC are regular (see Proposition 6.9). 

Emphasize that methods used in [9, 3] essentially use condition Q G T^[0,1] (g) (i.e. the 

fact that the Fourier coefficients belong to l‘^{TLj) and most likely could not be applied even 
to Dirac operators with L^-potentials Q. Note also that traditional methods of perturbations 
theory are also not applicable here since as opposed to the L^-case, the multiplication operator 
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by Q(s 1] 0 is neither B - compact nor even subordinated to the unperturbed 

operator 

The Riesz basis property for abstract operators is investigated in numerous papers. Due to 
the lack of space we only mention [19, 36, 37, 1], the recent paper [47] and the references therein. 

The main results of the paper including Theorem 1.1 were announced in [28] (partially with 
proofs). After appearance of [28] there appeared the paper by A.M. Savchuk and A.A. Shka- 
likov [44] where Theorem 1.1 was proved for the 2x2 Dirac operator. Note that approaches 
in [28] and [44] substantially differ. Moreover the case of Dirac type operators (6i + 62 7^ 0) has 
interesting features (see e.g. Corollary 1.2) and turns out to be more complicated. 

The paper is organized as follows. In Section 2 we prove the existence of triangular trans¬ 
formation operators for equation (1.8)-(1.9). In Section 3 we apply these operators to obtain 
asymptotic formulas for solutions to equation (1.8). In turn, these formulas are applied in 
Section 4 to obtain the following asymptotic formula 

An = -|-0(1), as re —>■ 00, re G Z. (1-11) 

for the eigenvalues A = {AnjnGZ of the operator Lc,d{Q) with regular BC. In Section 5 we 
present certain necessary and sufficient algebraic conditions for equations (1.10) to determine 
strictly regular BC. In particular, we show in Proposition 5.5 that if a := —bi/b 2 0 Q and a = 0, 
5c, d G M \ {0}, then BC (1.10) are strictly regular if and only if 

d 7 ^ —(a-|-1) (15c|a““) “+D ( 1 - 12 ) 

So, under the above restrictions condition (1.12) gives the algebraic criterion of the strict reg¬ 
ularity of boundary conditions (1.10). In Section 6 we prove our main results on Riesz basis 
property of the root vectors system of the operator Lc,d{Q) (Theorem 1.1 and Proposition 6.9). 
Finally, in Section 7 we apply Theorem 1.1 to prove the Riesz basis property with parentheses 
for the dynamic generator of the Timoshenko beam model (see e.g. [53, 22, 49, 58, 57, 56]). 

Notation. Let T be a closed operator in a Hilbert space Sj. Denote by p{T) the set of regular 
points of T; cr{T) = C \ p{T) and crp(T) denote the spectrum of T and the point spectrum of T, 
respectively. 

For the eigenvalue Aq G o'p{T) denote by ma{Xo) and mg{\o) the algebraic and geometric 
multiplicities of Aq, respectively. Recall that mg{Xo) = dim(ker(L — Aq)) and ma(Ao) is a 
dimension of the root subspace corresponding to Aq. 

Or(2:) C C denotes the disc of radius r with a center 2;. 

(•, •) denotes the inner product in C"'; denotes the set of re x re matrices with complex 

entries; /n(G C”^”) denotes the identity matrix. 


2. Triangular transformation operators 


2.1. The Banach spaces Xi and Xoo. Following [32] denote by Xi := Ai(D) and X^o ■= 
Xoo(D) the linear spaces composed of (equivalent classes of) measurable functions defined on 
D = {(x, t) : 0 < t < X < 1} satisfying 


11/llxi := esssup / \f{x,t)\dx < 00 , 
4G[0,1] Jt 


Yoo := esssup / \f{x,t)\dt < 00, 
xG[o,i] Jo 


( 2 . 1 ) 


( 2 . 2 ) 
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respectively. It can easily be shown that the spaces Xi and equipped with the norms (2.1) 
and (2.2) form Banach spaces that are not separable. Denote by Xi^q and X^o^ the subspaces 
of Xi and X^o, respectively, obtained by taking the closure of continuous functions / G C(D). 
Clearly, the set C'^(D) of smooth functions is also dense in both spaces Xi^o and Xoo,o- 
To motivate appearance of the spaces Xi and X^o consider a Volterra type operator 

N: r N{x,t)f{t)dt (2.3) 

Jo 

with a measurable kernel and denote by ||A^||p := ||.^||i,p[o,i]->-i,p[o,i] the Lp-norm of the 

operator N provided that it is bounded. The following simple lemma (cf. [32]) sheds light on 
appearance of the spaces Xi and Xoo- 

Recall that a Volterra operator in a Banach space is a compact operator with zero spectrum. 

Lemma 2.1. LetN{-,-) G Vi(D)nVoo(D) and generate the Volterra type operator (2.3). Then: 

(i) The operator N is bounded in LP[0, 1] for each p € [1, oo] and 

ll"ll.<IAIIx'f,n,-IIA'llVm)' (2-4) 

Moreover, 

IliVlIi = l|iV|Ui(0), l|A'lloo = ||iV||A^. (2.5) 

(ii) If N{-, ■) G Vi^o(^)nXoo,o(^2), then N is a Volterra operator in IT[Q,l] for eachp ^ [IjOoj- 

Proof, (i) The relations (2.5) are well known (see e.g. [33]) and can easily be proved. Combining 
the M. Riesz’s interpolation theorem with relations (2.5) yields 

IliVllp < IIA'Ilf” • jliVllL-'-'" = l|iV|lf% ■ IDIlVm. p e (2.6) 

which proves estimate (2.4). 

(ii) Since G Vi^o(^^) CX(X)^o(^^)) there exists a sequence A^fc(-,-) G (7^(11), /c G N, such 

that limfc_^oo(||^^ - A^fclUi(o) + 11^ - ^felUoo(n)) = 0. In accordance with (2.4) 

||iV - NkWp < [jiV - • IIV - p G [1, cx)]. (2.7) 

Since G (7^(12), the operator Nj^ of the form (2.3) is a Volterra operator. It follows 

from (2.7) that the operator N of the form (2.3) is the uniform limit in L^[0,1] of the Volterra 
operators and is therefore itself a Volterra operator. □ 

The following simple properties of the class Xoo,o(^2) will be useful in the sequel. 

Lemma 2.2. For each a G [0,1] the trace mapping 

ia ■■ (7(D) (7[0,a], ia{N{x,t)) := N{a,t), (2.8) 

originally defined on (7(D) admits a continuous extension (also denoted by ia) as a mapping 
^cx),o(^) L^[0,o] from Voo,o(^^) onto L^[0,a]. 

Proof. Let N{-, •) G Voo^o(^) and let Vfc(-, •) G (7(D) be a sequence approaching N in Xoo(D). 
It follows from definition (2.2) of the norm in Voo(D) that 

\Nk{a,t) - Njnia,t)\dt < \\Nk - NmWx,,^ ^ 0 as n,m^oo, (2.9) 

i.e. the sequence Nk{a,-) is a Cauchy sequence in L^[0,a]. Thus, there exists /„(•) G L^[0,a] 
such the ll/a — Nk{a, OllLbo.a] —>■ 0 as /c —>■ oo. We put N{a, •) := /a(-) and extend the mapping 
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ia to the space X^ofi by setting ia : X{-,-) —)■ /a(-) = N{a,-). It is easily seen that this 

extension is well defined. Indeed, if Nk{-,-) is another sequence approaching in X^o- 

Then limfc^oo \\Nk - Nk\\x^ =0 and 

lim \\Nk{a, •) - iVfc(a, OllLqo.a] < ,1™ ll^fc - ^k\\x^ = 0. (2.10) 

k^oo k^oo 

Hence lim^-^oo ll-^(O) 0 “ Xk{a, Olliqo.a] = 0 and the extension ia is well defined. □ 

Going over to the vector case we introduce the Banach spaces 
^2X2 _ ^ ^2X2 ;^^x2 _ X^^ (Q) := Xoo(I2) ® (2.11) 

consisting of 2 x 2 matrix functions / = {fjk)‘j k=i "'^itb entries from Xi and Xqo, respectively, 
and equipped with the norms 


II/IIyi := ||/||yi®C2x 2 := max{||/jfc|Ui ■ j,k e {1,2}}, (2.12) 

Yoc := II/IIyooOC2x2 := max{\\fjk\\x^ ■ j,k e {1,2}}. (2.13) 

We also put 

■= ■= ^i,o(b!) 0 C2x2 and := Xoo,o(fl) ® C2x2. 

Further, equip the space T^([0,1],C^) := L^[0,1] (8) of vector functions with the following 
norm 

ll/llp := ||col(/i,/ 2 )||p := ||/i||p+||/ 2 ||p, pG[l,oo]. (2.14) 

where ||/j||p := ||/j||lp[o,i], J e {1,2}. 

With each measurable kernel {^Njk{-, O)^fc=i oae associates a Volterra type operator 


N 


r N{x,t) 

Jo 


hit) 

hit) 


dt = 


/iVii(a:,t) Ni2ix,t)\ fhit)\ 

\N2iix,t) N22ix,t)) \f 2 it)) 


dt. 


(2.15) 


Let us set ||A^||p := ||lV||p^p[o,i](g)C2-)-Lp[o,i](g)C2) P € [l,oo] provided that the norm is bounded. 

Lemma 2.3. Let N{-,-) = (iy 7 fc(‘, O)^ fc=i) € W^x2(f2) nX^^(n) and generate the Volterra type 
operator by formula (2.15). Then: 

(i) The Volterra type operator N is a bounded operator in L^([0,1],C^) for each p G [l,oo] 
and 


p — 


Moreover, 


1 — 




iixr 

^00 

(2.16) 

and 

11^^1100 = ||iV|lxix2. 

(2.17) 


(ii) If G W^Q^(n) n X^Q(n), then N is a Volterra operator in LP([0, 1],C^) for each 

P G [ 1 , 00 ]. 


The proof is similar to that of Lemma 2.1 and is omitted. 

Next we demonstrate that the assumption N(-, ■) G Xi^o(12) G X^ofti^) in Lemma 2.1(ii) is 
essential for the operator At to be a Volterra operator. 
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Proposition 2.4. Let k{-) € L^[0,1], k{s)s ^ G L^[0,1], and let Aikia) := Jq k(s)s° ^ ds be 
the Mellin transform ofk[-), a a Cr ■= {z € C : Rez > 0}. Then the Volterra type operator 

K-. f^ - [\(-)fit)dt. (2.18) 

X Jo \xj 

has the following properties: 

(i) The operator 1C is bounded in ^^[0,1] for each p G [1, oo]. 

(a) Its point spectrum is given by ap{IC) = range{A4.k{-)) = {fo ds : a G Cr}. In 

particular, K, is not compact. 

(Hi) If k{-) > 0, then the spectral radius of K, is equal to its norm \\IC\\p in each L^[0,1]. 


Proof, (i) Let us check that N{x,t) = ^k{-) G Xi(r2) n Xoo{Cl). Indeed setting t/x = s one 
easily gets 


\N{x,t)\dx = J —\k^—j\dx = J^ s ^|A:(s)|ds< J s ^|A:(s)|ds, 


and 


[ \N(x,t)\dt = — [ Ifef —[ —|A;(s)|a;ds 
Jo X Jq \xJ Jo x 

= f |A:(s)|(is< f 

Jo Jo 

The boundedness of /C in L^’[0,1], p G [1, oo], is now implied by Lemma 2.1. 
(ii) Clearly, fa = G L^[0,1] for a G and 


{ICf, 


Mx) = - [ 
X Jo 


"" k 


1 

e-^dt=- / k{s){xs)^-^xds = Mk{a)x^-\ 

X Jo 


Due to the assumption k{-)fo{-) G L^[0,1] one has 

|T4(q;)| < [ \k{s)s°‘~^\ds < [ ds =; M|fc|(0). 

Jo Jo 


(2.19) 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 


Since the function Mk{-) is holomorphic and bounded in C^, it might have at most countable 
(discrete in Cr(—1)) set of zeros. For the rest of ath is the eigenvector of R in L^[0.1] 
belonging to a non-zero eigenvalue c(a), i.e. x" G ker(/C — c{a)). Hence /C is not a compact 
operator. □ 


2.2. Transformation operators. 

Theorem 2.5. Let Q12, Q21 G -^^[0, Ij- Assume that e±(-; A) are the solutions of the system (1.8) 
corresponding to the initial conditions e±(0;A) = (j_^)- Then e±(-;A) admits the following 
representation by means of the triangular transformation operator 


e±{x-, X) = {I + K^)e^{x-, X) = e^{x-, X) + I K^{x,t)e^{t-, X)dt, 


±1 


(2.23) 


where 


( pibiXx \ 

^))j,k=v ( 2 - 24 ) 
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and nX^Q(r2). In particular, the operator : f —)■ K^{x,t)f{t)dt, is 

a Volterra operator in each L^lO, 1], p G [l,oo], hence a{K^) = {0}. 

Our further considerations will substantially be relied on the following result which is a special 
case of [33, Theorem 1.2] where the general case of n x n system (2.26)~(2.28) with the matrix 
B = B* € and Q G C[0,1] <Si was treated. 

Proposition 2.6. [33] Let Q = codiag((5i2, (52i) G C^[0,1] Then the boundary value 

problem 

B-^D^K^{x, t) + DtK^{x, t)B-^ + iQ{x)K^{x, t) = 0, (2.25) 

K^{x,x)B~^ — B~^K^{x,x) = iQ{x), xG[0, 1], (2.26) 

{x,0)B Q J =0, X G [0,1]. (2.27) 

has the unique solution •))^fc=i ^ C'^(n) Moreover, is the 

matrix kernel of the transformation operator (2.23). 

The proof of this result in [33] is divided in two steps. At first it is proved solvability (and 
uniqueness) of the certain auxiliary boundary value problem which in conformity to the 2 x 2-case 
reads as follows 

B~^DxR{x, t) + DtR{x, t)B~^ + iQ{x)R{x, t) = 0, (2.28) 

R{x,x)B~^ — B~^R{x,x) = iQ{x), xG[0, 1], (2.29) 

i?ii(x,0) = i?22(x,0) = 0, xG[0,l]. (2.30) 

where R{x,t) = {Rjk{x,t)^‘^. Let us recall the corresponding statement from [33]. 

Proposition 2.7. Let Q = codiag(Qi 2 , Q 21 ) G (^^[0,1] (8> Then the auxiliary prob¬ 
lem (2.28)~(2.30) has a solution R G (7^(12) <8> Moreover, it is unique in A^^(f2). 


Our first auxiliary result reads as follows. 

Proposition 2.8. Assume thatQ,Q G C^[0, 1](8>C^^^ and ||<5||li[o,i]®C2x2) IIQIIli[o,i] ^ x. Then 
there exists a constant C = (^(r, 61,62) such that 

ll-R - ~ QllLi[o,i](giC2x2. (2-31) 


Proof. We put 


and 


Qj := b- and ;= — = -p-, 

a? bk 


j,ke { 1 , 2 }. 


ijk{x,t) = 


{okX - ajt){ak - Oj) \ j ^ k, 


X — t^ j = k. 

Let us rewrite boundary value problem (2.28) - (2.30) in the scalar form 

{DxRkk (^; ^) “ 1 “ DfRkk (^; 0 ) — iQkj ^) ; ^ ^ ; 2 }, 

^jDxRjki^'it) Clf^DtRjki,^')^) — '^Qjk{^)Rkk{^i^)^ k G {1,2}, j ^ fe, 

Rjk{x,x) = xG[0,1], 

dk dj 

i?ii(x,0)=0, R22{x,Q)=D, xG[0, 1]. 


(2.32) 

(2.33) 

(2.34) 

(2.35) 

(2.36) 

(2.37) 
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The system (2.34) - (2.35) is hyperbolic. Integrating the Goursat problem (2.34) - (2.37) along 
characteristics we arrive at the following equivalent system of integral equations 


Rkk{x,t) = - - Qkji^Rjki^,^ - X+ t)d^, k € {1,2}, 


dn 


(2.38) 


J Jx-t 


Rjk{x,t) = - Qjk' 


dj 


dk dj 


— f Qjk{C)Rkk{C,Kjk{C - x)+ t)d^. (2.39) 

«A; 


Rere j,k £ {l,2},^j ^ k 

The functions Rjk satisfy the same system (2.38)-(2.39) with Qjk instead of Qjk, j, k € {1, 2}. 
Next we put 


Rjk Rjk Rjk: Qjk Qjk Qjki ji ^ ^ {Ij 2}. 


and 


^ px ^ px 

5 jfc(x) := / \Rjk{x,t)\dt, 3 jkix) := \Rjkix,t)\dt, j,k£ { 1 , 2 }. 

Jo Jo 

Making use the change of variables ^ = u, ^ — x + t = v, we obtain from ( 2 . 38 ) that 

^ px ^ px px ^ 

5 ii(x) := / \Rii{x,t)\dt < \bi\ / dt \Qi 2 { 0 R 2 i{C,C - x + t)\dC 
Jo Jo J x—t 

rx rx 

+ \bi\l dt \Qi 2 {OR 2 i{^,^ - X+ t)\dC 
Jo J x—t 

px ^ pu px ^ pu ^ 

<\bi\ / 1^12(14)1^11 / 17221(11,14)1^1; + 1611 / |Qi 2 (i 4 )|di 4 / |.R 2 i(i 4 ,i’)|di’ 

Jo Jo Jo Jo 

px ^ px ^ ^ 

= \bi\ I \Qi2{u)\J2i{u)du + \bi\ I \Qi2{u)p2i{v) du. 

Jo Jo 

Similarly, it follows from ( 2 . 39 ) and ( 2 . 41 ) 

52 i(x) = f \R2i{x,t)\dt < - - - - 7 f Q2 i(^^^ 

.10 I«1 — 0^21 Jo V 

px pa 

H-/ dt 

02 Jo Je .2 




CLl — 02 

1^21(0-^21 (?, 4121 (C - X) + t)|(i^ 


+ ■ 


"I px px 

- dt \Q 2 i{C)R 2 i{^,i^ 2 i{^ - x) +t)\dC. 

12 JO Ji 2 i{x,t) 


(2.40) 


(2.41) 


(2.42) 


(2.43) 
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Making use the change of variables i = u, — x) +1 = v, we obtain 

^ 2 i{x) = — \Q 2 i{u)\du-\ -/ \Q 2 i{u)\du \Rii{u,v)\di 

02 J 02 J J^(u-x) 

ai—0.2 ai—^2 ^2^ ^ 

1 


01-02 
cx 


+- 


02 


/ X ^ ru ^ 

\Q 2 i{u)\du / \Rii{u,v)\dv 

02 ^ ^ 


^ j^X ^ ^X ^XL 

< — / |Q2 i(^‘)|c^'0 H-/ |Q2 i('o)Mo / 

02 Jo 02 Jo Jo 

-j j* ^ n 

-/ |Q2 i(w)Mo / 

02 Jo Jo 

-j ^x ^ ^X ^ f*X 

— / \Q2i{t)\dt^ -/ Jii(rt)|Q2i(o)|du H-/ 5ii(u)|Q2i(ii)|d?.i- 

02 Jo 02 Jo 02 Jo 


Let 


C'j — ||flji||L°°[o,i] — l|.Rii||xoo(D)) J € { 1 ) 2 }. 

Estimate (2.42) with account of (2.45) yields 

5 ii(x) < |6 i|C'2||Qi2||li + \bi\ [ ^ 2 i{u)\Qi 2 iu)\du. 

Jo 

Combining this inequality with (2.44) implies 

( 1 C \ 1/*^^ 

-^- -] |Q2i(o)|d'oH-/ fJii('o)|( 52 i('o)|do 

02 02 / Jo 02 Jo 

cx PU 

< C'(||Q2i||i,i + |ftl^2|C'2||<3l2||Ll ■ ||Q2 i||l 1 + 1^1^21 / |<32l(o)|(it( / 3i21 (t)|(5l2(0 

Jo Jo 


(2.44) 

(2.45) 

(2.46) 


<6*3 


2l||l,l + 


px px 

IIQi2||i,i) + |fel&2| y ^ 2 lit)\Ql 2 {t)\dt \Q 2 l{u)\du 


<Cs 


2 i||l 1 + II<3i2||l1 ) + 


\bib2\ ■ \\Q2 i\\l^ [ '^2i{t)\Qi2it)\dt, (2.47) 

Jo 


where 


C{:= 


l + Ci 


02 


and Cs := max{C(, |6 i 62|C'2||Q2 i||li}- 
Applying Cronwall’s lemma to this inequality implies 

52 i(x) < Cs f||(g2i||Li + ||Qi 2 ||lO exp 


|6i62|-||Q2i||li \ Q12mty 
Inserting this inequality in (2.46) we arrive at the inequality 
3111(0;) < |6i| f||Ql2||Ll + IIQ21IIL1) (<32 + (3311(51211^1 exp (I6162I • ||(5 i 2 ||l 1 ' ||(52l||l,l 


(2.48) 


(2.49) 

Similar reasoning leads to similar estimates for fJi2 and 322- Combining these estimates 
with (2.48) and (2.49) we arrive at (2.31). □ 

Proposition 2.9. Let Q = codiag((5i2) (52i) S A^[0,1] (g) Then the system of integral 

equations (2.38) - (2.39) has a unique solution R = {Rjk)‘j k=i belonging to A^Q(fl). Moreover, 
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Further, let Qn = codiag((5i2,n, Q2i,n) £ C'^[0,1] 0 he any sequence approaching Q in 

-L^[0, l]-norm and let Rn = {Rjk,n)‘j k=i ^ be the corresponding system of solutions 

of the problem (2.28)-(2.30) with Qn instead of Q. Then there exists a constant C = C{Q, bi, 62) 
not depending on n and such that the following estimates hold 


“1“ ^jk,n\\Xoo — ^IIQ ^ ^ ^ ^ ^ (2.50) 

Proof, (i) Choose sequences {(5i2,n}nGN) {Q2i,n}neN C (7^)0,1] such that 

IIQ12 ~ Ql2,n||Ll + IIQ2I ~ Q21,n\\L^ —>■ 0 aS U —>■ OO. (2-51) 

By Proposition 2.7, for each pair Qn = {Qi2,nj Q2i,n} there exists the unique matrix solution 
Rn = {Rjk,n)‘j k=i ^ C'^(n) (g) n e N of the system (2.38)-(2.39). It follows from (2.51) 

and (2.31) that there exists R = {RjQ'j S X^^(n) such that 

lim \\Rjk,n ~ ^jkWxoad^) ~ ^ {1)^}’ (2.52) 


Let us show that {Rjki'i ■)}‘j k=i satisfies the system (2.38)-(2.39). Let for instance equation 
j 7^ k. Writing down equation (2.39) for Rjk,n{'^') and integrating it with respect to t from 0 
to X one gets (cf. (2.44)) 


j-x i i /■* /■“ 

/ ^jk,n{^R)dt — / Qjk^n{u)du / Qjk^niu)du I Rkk,n{^U,v') 

Jo a? J a? J J^(u-x) 

"■j 


dv. 


(2.53) 

It follows from estimate (2.48) that Wuin^oo ^kk,n{u-,v)dv = Rkk{u,v)dv for any pair 
vi,V 2 € [0,1]. Therefore and due to (2.52) the dominated convergence theorem applies as 
n —>■ OO in (2.53) and gives 


j-x i i 

/ Rjkix,t)dt = -/ Qjk{u)du -/ Qjkiu)du Rkk{u,v) 

Jo J J—J ^^{u—x) 


dv 


a^-aj 

[ — - — Qjki ——/ Qjk{ORkk[f,,i^jk{i-x) + t)df 

Jo afc - “i \ ak- aj ) Uk J^.Xx,t) 


ai^-aj 




dt. (2.54) 


The latter is equivalent to (2.39). The equations for Rjj(-,-), j € {1,2}, is obtained similarly. 

(ii) Since the sequence Qn{-) approaches Q{-) in L^-norm, it is bounded, ||Qn||Li(giC2x2 < Ci = 
Ci{Q,B), n € N. Therefore Proposition 2.8 applies and gives 

[[.l^jA; dijk,n\\Xao — C*IIQ Qn\\L^^C^^^t j) k G {1, 2}. (2.55) 

Next we prove similar estimate in Xi(ll)-norm. We let 

Bjkn — dJ-jk dJjkm Qjkn — Qjk Qjkm j: k G {1, 2}, n G N. (2.56) 

First we prove estimate (2.50) for the case j ^ k. To this end we note that 

[ Qjkn(— - —^dx=[ Qjkn{ijk{x,t)) dx = ——^ / Qjkn{u)du. (2.57) 

Jt \ afc o.j J ak Jt 

where f,jk{-, •) is given by (2.32). 
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Further, note that Rjkn{-, •) satisfies equation (2.39) with Qjkn in place of Qjk, j,k € {1,2}. 
Taking difference of this equation (2.39), then integrating the difference with respect to x € [t, 1], 
and making use the change of variables u = v = — x)Kik + t = — x)^ +1, we obtain 




It 


\Rjk,n{^,t)\dx = / \Rjkix,t) - Rjk,n{x,t)\dx 


< 


\ak\ 


I ^ 1 rijkO -,!) /■('"-*)^+1 

/ \Q {u)\du + — 

Jt \aj\ Jt 

1 fijkiu) _ 


\ak\ 


I \Qjk,n{'ll')\ dlt- r / dx I \Qjk(.C^Rkk,n(.^^ I^jki.^ x') t')\d^ 

Jt \(^j\ Jt Jijk(x,t) 

T"! r / dx I \Qjk,n{^')Rkk,n(.^il^jk{^ x) t')\d^ 

/■(’'-*) ^+1 

I \QjkinJ)Rkkjnilli l^^\du 

J V 

1 rijk(k,t) |■{v-t)^+l ^ 

7—7 dv |<5jfc,n(«)^fcfe,n(«,I^)M« 

\^j\ Jt j V 

7^ / \Qjk,n{u)\du +J— I dv I \Qjk{u)Rkk,n{'^,v)\du 

\^k\ Jt Jt Jv 

\^j \ Jt JV 

\Qjk,n(u)\du +\Qjk{u)\du \Rkk,n{u,v)\dv 

\Qjk,n{'^)\du / \Rkk,n{u,v)\dv 

J V 


< 


1 

|nfc| Jt 


— I^fcl ■ IIQjA:,n||Ll[t,l] T \bj\ ' IIQjfc||Ll[t,l] ' ||-^fcA:,n||Aoo 

+ \bj\ ■ ||(5jfc,n||l,i[t,l] ■ ||.RA:fc,n||Xoo ■ (2.58) 


Here we use simple inequalities t,jk{^,t) < 1 and {v — t)— + 1 < 1. The latter holds since t > v 
and ajUk < 0. It follows from (2.58) with account of definition (2.1)-(2.2) that 


^ \bk\ ' l|l?jfc,n||Ll[0,l] 

l^il ^IIQjfclll,! ■ ||l^fcfc,n||Aoo T I|l5jfc,n||l,i ■ ||.Rfcfc,n||Xoo^ , j 7^ k. (2.59) 

On the other hand, estimate (2.55) implies lim„_).oo \\Rjk — Rjk,n\\xoo — 0- Therefore there exists 
C 2 > 0 such that imLK{\\Rjk^n\\xoo ■ J^k G {l,2},n G N} < C2. Combining this estimate with 
(2.55) yields the following estimate 

ll-^jfc Rjk,n\\Xi ^C* 3 ||Q Qyi||£l(g|f; 2 x 2 , J, A:G{1,2}, (2.60) 

with a ceratin positive constant 6*3 > 0 not depending on n G N. Combining this estimate 
with (2.55) implies (2.50) for j 7^ k. 

(hi) Going over to the case j = k we start with equation (2.38) and similar equation for 
Rkk,n{-, •) which holds with Qjk,n in place of Qjk, j,k G {1,2}. Taking difference of this equation 
and (2.38), then integrating the difference with respect to x G [t, 1], and then making use the 


14 


ANTON A. LUNYOV AND MARK M. MALAMUD 


change of variables ^ = u,^ — x + t = v, obtain as 

/■i 


< 


1*1 nx ^ 

I dx I 

Jt Jx—t 


i: 


aj\ I \Rjj^n{x,t)\dx = \aj\ j \Rjj{x,t) — Rjj^n{x,t)\dx 


^-x + t)\d^+ dx \Qjk,n{ORkj,n{^,C - X+ t)\d^ 


'x—t 


rt ^ pi pv-i-i-i ^ 

= dv \Qjk{u)Rkj,n{'^,v)\du + dv \Qjk,n{u)Rkj,n{u,v)\du 

Jo Jv Jo Jv 

< dv \Qjk{u)Rkj,n{u,v)\du+ dv \Qjk,n{u)Rkj,n{u,v)\du 

J 0 J V J 0 J V 

pi pu pi pu 

= \Qjkiu)\du \Rkj,niu,v)\dv+ \Qjk,ni'^)\du \Rkj,niu,v)\dv. (2.61) 
Jo Jo Jo Jo 

It follows with account of definition (2.1) - (2.2) that 


fV — t+l 


W^jj dijj^nWxi ^ l^jl ^IIQjfcllLl ■ l|-^A:j,n||Aoo “t” IIQjA:,n||Ll ' l|-^A:j,n||Aoo^ J J ^ ( 2 . 62 ) 

Since \\Rjk,n\\x^ < C 2 for n € N, this estimate together with (2.55) leads to the estimate (2.50) 
with j = k. □ 


Lemma 2.10. Let Qi 2 ,Q 2 i G -^^[0,1] and letR{-,-) = {Rjk{■,•))] k=i ® solution of the system 

of integral equations (2.38)-(2.39). Then R{-,-) £ X^Q^(fl) nX^^Q(H) and the operator 


R : 





r (R\x{x,t) Ri2{x,t)\ (fi{t)\ 
Jo \R2lix,t) R22ix,t)) \f2{t)J 


(2.63) 


is a Volterra operator in LP[0, 1] (g) for each p £ [1, 00]. 


Proof. Let Qn = codiag(Qi2,rn (52i,n) £ 1] ® be any sequence approaching Q in 

L^[0, l]-norm and let ii„ = {Rjk,n)‘j k=i ^ C'^(^^) n G N, be the corresponding system of 

solutions of the problem (2.28)-(2.30) with Qn instead of Q. 

Since Rn = {Rjk,n)'j k=i ^ smooth kernel and hmn_>.oo ||Q — Qn||Li[o,i]®C2x2 = 0, it follows 
from Proposition 2.9 (see estimates (2.50)) that R{-,-) ^ Therefore by 

Lemma 2.3, i? is a Volterra operator in LP[0, 1] g) C^, p £ [1, 00]. This completes the proof. □ 

Proof of Theorem 1. Let P^ = diag(Pj^,P^) be a diagonal matrix function with entries Pp £ 
L^[0,1], j £ {1,2}. Define the convolution operator 

P^ : f^ [ P^{x -t)f{t)dt, / = col(/i,/2) G L^[0,1] gC^. (2.64) 

Jo 

Let R{x,t) = {Rjkix,t))‘j be the solution of the system of integral equations (2.38)~(2.39). 
Starting from the operator I + R and following the reasoning of [33, Theorem 1.2] we define the 
operator by the equality 

I + K^ = {I + R){I + P^). (2.65) 

The latter means that the kernel K^{-, •) of is given by 

/ X 

R{x, s)P^{s — t)ds. 


( 2 . 66 ) 
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Let US show that K^{-, •) is the kernel of the transformation operator, i.e. representation (2.23)- 
(2.24) holds. First, we choose P^[-) so that ■) will satisfy condition (2.27), i.e. 

aiKflix, 0) + a2Kf^{x, 0) = 0, j € {1, 2}. 

Inserting representation (2.66) for K^{-, •) in these relations leads to the following system of 
Volterra type integral equations 


aiP^{x) + J^[aiRu{x,t)P^{t) ± a 2 Ri 2 {x,t)P 2 {t)]dt = ^a2Ri2{x,0) =: gfix), 
a 2 P^{x) + J^[aiR 2 iix,t)P^{t) ± a 2 R 22 ix,t)P^{t)]dt = -aiR 2 i{x, 0 ) =: gf{x). 


(2.67) 


Here the relations (2.37) have been taken into account. It follows from equation (2.39) that the 
functions Rjk{x,0), j,k G {1,2}, are well defined. Moreover, the estimate (2.58) ensures that 
9j{') £ 1], J £ {1, 2}. 

On the other hand, by Lemma 2.10 the operator R of the form (2.63) is a Volterra operator in 
L^[0,1]. Therefore system (2.67) is the system of Volterra equations in L^[0,1] with respect 
to col{aiPf*^(-),a2P^(-)}, hence has the unique solution col{aiPf*^(-),a2H2^(')} ^ 1] 

Further, choose a sequence Qn = codiag((5i2,n, Q2i,n) £ 1] (8> approaching Q in 

L^[0, l]-norm. Then according to Proposition 2.7 there exists the corresponding sequence of 
matrix solutions Rn = {Rjk,n)‘j k=i ^ C^(f2) (S> of the problem (2.28)-(2.30) with Qn 
instead of Q. Moreover, by Proposition 2.9, the estimate (2.50) holds, hence Rn approaches R 
in Xi and Xqo norms. Choose a sequence P^ = diag(P^^, P^^) of diagonal matrix functions 
with entries Pfn(-) G C^[0,1], J G {1,2}, re G N, and assume that P^(-) satisfies the following 
system of Volterra integral equations 


aiP^^(x) + fg[aiRn,n(x,t)P{‘^n(^) ± a2Rl2,n(x,t)P:^n(^)]^^ = =FCi2Rl2,n(x,0) 

^0^2.P2,ni^) + /o'[“l-^21,n(a;,l)L’i|ii(l) ± a2R22,n{x,t)P2^nW\dt = -aii?21,n(a:, 0) =: 5^„(x) 


Next we define the kernels K^{-, •) by setting (cf. formula (2.66)) 


Kn{x,t) = Rn{x,t) + ^nix - t) + / Rn{x, s)Pn {s - t)ds, 


re G N. 


( 2 . 68 ) 


(2.69) 


Clearly, K^{-, •) G C'^(n) and in accordance with [33, Theorem 1.2], it is the unique solution 
of the boundary value problem (2.25)-(2.27). Note for instance, that condition (2.26) for the 
kernel Kn{-,-) is satisfied since Rn{-,-) satisfies this condition, Kn{x,x) = Rn{x,x) + P,^(0), 
and the matrix P^(0) is diagonal. 

Further, by Proposition 2.6, iL)}(-,-) is the kernel of transformation operator for equa¬ 
tion (4.1)-(1.9) with Qn in place of Q, i.e. the solution e±^„(-;A) of this equation satisfying 
the initial condition e±_„(0; A) = (j_^) admits a representation 

e-±,n{x',\) = {I + Kn)e\.{x-,\) = e\{x-^\) + [ Kn{x,t)e^{t;X)dt, re G N. (2.70) 

Jo 

Our aim is to pass to the limit in (2.69) and (2.70) as re ^ oo. It follows from (2.58) with 
t = 0 and the estimate (2.50) that 


lim ( 

n—)-oo 


tn - IIl 1[0,1] + htn “ 11^1 [0,1]) = 0- 
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Combining this relation with Proposition 2.9 we obtain from (2.67) and (2.68) that 


.^2PU-) 


,± / X ) ~ + Pn] 


1") 


. 92 , J 

\92 


Ol-Pf (•) 

a2PH-) 


(2.71) 


Further, setting := Pn — P we derive from (2.66) and (2.69) 

f \K^{x,t)\dx = [ \K^{x,t) — K'^{xp)\dx < f \Rn{x,t)\dx+ f \P^{x — t)\dx 

Jt Jt Jt Jt 

+ [ —^)Ms f \Rn{x, s)\dx + f \P^{s — t)\ds f \R{x,s)\dx 

J t J S J t J S 

< Il-Rn||xi(r2) (l + II-P^IIli[o,i]) + II-P^IIli[o,i] (l + ll-R||xi(o)) • (2.72) 

On the other hand, by Proposition 2.9, lim„_>.oo = hm„_>.oo ||7?n — R\\x‘^^'^ = 0, and 

due to (2.71) lim„_>.oo ||P’^||li[o,i] = 0- Combining these relations with (2.72) yields 

lim IliFi^ll v 2 x 2 = lim IliF^ — 771(^11 ^ 2 x 2 = 0. (2.73) 

n^oo n^oo 

The latter means that G X^q^(O). In just the same way one proves the relation 

lim ||iFi)^|| ^ 2 x 2 = lim 1177^ — ^ 2 x 2 = 0. (2.74) 

n-)-oo " "■ ""^ 0 ° n^oc " "■ ^ ^ 


Using relation (2.74) we can pass to the limit as re —>• 00 in formula (2.70) and arrive at the 
required formula (2.23). □ 

Remark 2.11. (i) For Dirac 2x2 system {B = diag(—1,1)) with continuous Q the triangular 
transformation operators have been constructed in [26, Ch.10.3] and [35, Ch.1.2]. For Q G 
7y^[0,1] ( 8 ) it is proved in [2] by an appropriate generalization of the Marchenko method. 

(a) Let J : f ^ Jq f{t)dt be a Volterra operator on L^[0,1]. Note that the similarity of 
Volterra operators given by (2.15) to the simplest Volterra operators of the form B ® J acting 
in the spaces L^[0,1] (8> has been investigated in [33, 43]. 


3. Asymptotic behavior of solutions 

v 2 


Let K^{x,t) = (iFi^(a:, t))^. be the kernel of a triangular transformation operator con¬ 
structed in Theorem 2.5 (see formulas (2.23)-(2.24)). To state the next result we put 


R%--=‘^~\K7k + Kik)^ J, fee {1,2}, 


jk 


jk) 


(3.1) 


and let 


SmO <!>.(■.>)), 4.(0, A) = /„ (3,2) 

be a fundamental matrix solution of the system (1.8). Here 4>fc(-, A) is the fcth column of 4>(-, A). 
Our investigating of the perturbation determinant relies on the following result. 
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Proposition 3.1. Let Q € L^[0,1] and let j,k € {1,2}, be the entries of the 

fundamental matrix solution (3.2). Then the functions admit the following representa¬ 

tions 


Jo Jo 


ipi2{x,X) = 


rx rx 

/ Rfi{x,t)e''^'^^^dt + / Rf2{x,t)e^^^^^dt, 

Jo Jo 

px px 

T2i{x,\)= / + / Rf2ix,t)e^^^^^dt, 

Jo Jo 

px px 

T22{x,X) = + / Rfi{x,t)e''^~^^^dt + / R22{x,t)e^^^^^dt, 

Jo Jo 


where n Xoo,o(ll), J, k G (1, 2}. 


(3.3) 

(3.4) 

(3.5) 

(3.6) 


Proof. Comparing initial conditions and applying the Cauchy uniqueness theorem one easily 

gets 4>i(-,A) = 'a)) ~ + 6 -(S"^)- Inserting in place of e+(-;A) and e_(-;A) their 

expressions from (2.23) one arrives at (3.3) and (3.5). Relations (3.4) and (3.6) are proved 
similarly. □ 


Lemma 3.2. Let N{-,-) € Xoo^o(12), b € M. \ {0} and h > 0. Then the following asymptotic 
holds uniformly m x G [0,1] 


A^(x, —>■ 0 as A —)• oo, |ImA| ^/i. (3-7) 

Proof. By the definition of the space 2foo,o(12), the inclusion N G Xao,o{^) ensures that for any 
e > 0 there exists G C'^(n) such that 

px 



||1V - = esssup / \N{x,t) - N^{x,t)\dt < e. 

xe[o,i] Jo 


In particular, we get the following uniform estimate 


f 


{N{x, t) — N^{x, t)) e^’^^^dt 


^ ee 


|b|h 


X G [0,1], |Im A| ^ h. 


(3.8) 


(3.9) 


Since G C^(n), integrating by parts the integral Jq N^{x,t)e^^^^dt we obtain the following 
estimate uniformly in x G [0,1] with some C > 0 


f 


N^{x, t)e^^^*dt 


< j^, A 7 ^ 0, |ImA| ^ h. 


The desired formula (3.7) now directly follows from estimates (3.9) and (3.10). 


(3.10) 

□ 


Remark 3.3. ITe demonstrate that the assumption N{-,-) G Xoo^o(ll) 'Is important for the 
validity of the statement of Lemma 3.2. More precisely, we show that for certain N(-, ■) G 
Xcx 3 (fl) \ Xoo,o(ll) the pointwice convergence in (3.7) holds but is not uniform in x ^ [0,1]. 

Let N{x,t) = yfc(^) where k{-) satisfies the conditions of Proposition 2.4. Then 



ds — y 0 as A —y oo, A G R/^, (3.11) 
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for each x € (0,1] and A G II/j. However this convergence is not uniform in x G [0,1]. Indeed, 
since k{-) ^ 0, its Fourier transform k{-) does not vanish identically, i.e. there exists a G M 
such that k{a) ^ 0. Therefore for A G M big enough and x = a/X G (0,1) the right hand side 
of (3.11) is k{a) 7 ^ 0. 

In the sequel we need the following result on the asymptotic behavior of solutions of the 
system (4.1) in the strip 

:= {A G C : |Im A| ^ h}. 

Proposition 3.4. Let Q G L^[0,1] 0 Then for any h > 0 the following asymptotic 

relations take place uniformly in x G [0,1] 

ifjkix, X) = + o{l) as A —>■ 00 , A G H/j, A:G{1,2}. (3-12) 

Proof. The proof immediately follows by combining Proposition 3.1 with Lemma 3.2. □ 


Applying the same approximation procedure as has just been used in the proof of Lemma 3.2 
to the space L^[ 0 , 1 ] we obtain the following simple statement useful in the sequel. 

Lemma 3.5. Let g G L^[0,1] and c G C \ {0}. Then for any e > 0 there exists M = Mg > 0 
such that 

< e(eR«pA)+ |A|>M. (3.13) 

Further, consider the adjoint system 

-iB~^y'+ Q*{x)y = Xy, a;G[0,1], (3.14) 

and introduce its fundamental matrix solution 

= (£m) = 

Here Tfc(-, A) is the fcth column of '!'(•, A). Clearly, Proposition 3.4 holds for the matrix solution 
'L(-,A) as well. Hence (3.12) and similar relations for ipjki'^X) imply the following result. 



Corollary 3.6. Let /i > 0. Then for X — )■ 00 , A G the following asymptotic relations hold 

($j(-,A),Tfc(-,A)) =(5jfc + o(l), j,A:G{l,2}, (3.16) 

(ci>i(.,A),4>2(-,A)) = o(l). (3.17) 

Moreover, there exist constants M > 0 and Ci,C 2 > 0, such that 

0<Ci < |(^j(-,A),$,(-,A))| <C2 , AgH;,, |A|>M, jG{l,2}. (3.18) 

Proof. First let us evaluate (<l>i(-, A), $!(•, A)). Setting f{x) := a; G ffi, and noting that 

giaAa:^ a G M, is bounded for (x. A) G [0,1] x H/^, one easily deduces from the uniform asymptotic 

relations (3.12) 

||4>i(-, A)|p = (4>i(-, A), 4>i(-, A)) = j [^pii{x, X)(pii{x, X) + ipi 2 {x, X)(pi 2 {x, A)^ dx 


g*H(A-A)x ^ 


-26iImA _ 1 

dx =-—— -h 0 ( 1 ) 

-26iImA ^ ^ 


=/(—26ilm A) + 0 ( 1 ) as A ^ 00 , A G H/, 


( 3 . 19 ) 
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Clearly, there exists Ci,C 2 >0 such that 

Cl ^ |/(x)| ^ C 2 , 


l^l ^ 2\bi\h. 


(3.20) 


Combining (3.19) with (3.20) one proves (3.18) for J = 1. Relation (3.19) for j = 2 as well as 
relations (3.16), (3.17) are proved similarly. □ 


4. Regular boundary conditions 
Here we consider 2 x 2-Dirac type equation (4.1), 

- + Q(x)y = Ay, y = col(yi,y 2 ), xe[0,1], 

subject to the following general boundary conditions 

[/j(y) := ajiyi(O) +aj2]/2(0) + aj3yi(l) + aj4y2(l) =0, J € {1,2}. 


(4.1) 


(4.2) 


Denote by L := L(Q, Hi, U 2 ) the operator associated in R^([0,1];C^) with the BVP (4.1)-(4.2). 
It is defined as the restriction of the maximal operator Hmax = Cmax(Q) (1-5) to the domain 


dom(L) = dom(L((5, Hi, U2)) = {y € dom(Hmax) : Ui(y) = U2(y) = 0 }. 


(4.3) 


The eigenvalues of the problem (4.1)-(4.2) are the roots of the characteristic equation A(A) := 
det H(A) = 0, where 


- ^Ci( 4 >i(.,A)) Hi( 4 . 2 (-,A))\ /uii(A) ni2(A) 

- Vc/2(^i(-,A)) H2(4>2(-,A))yl - V^2 i(A) n22(A) 


(4.4) 


Putting Ajk = j ^ j,^ _ g {1,...,4}, we obtain the following 

\CL 2 j 0 , 2 k J 

expression for the characteristic determinant 

A(A) = J 12 + + J32¥?1i(A) + Jl3(/9l2(A) + Ja 2 ^ 2 i{A) + Jia^ 22 {A), (4.5) 

where ^jk{A) := ^pjk(^,A). If Q = 0 then ifi 2 (x,X) = if 2 i(x,A) = 0 and the characteristic 
determinant Ao(-) becomes 


Ao(A) = J12 + ^ j^^gi62A_ 

In the case of Dirac system (B = diag(—1,1)) this formula is simplified to 

Ao(A) = J12 + J34 + J 32 e + Ji 4 e*^. 


(4.6) 


(4.7) 


Substituting formulas (3.3)-(3.6) at a: = 1 to (4.5) and taking into account (4.6), we get the 
following expression for the characteristic determinant. 

Lemma 4.1. The characteristic determinant A(-) of the problem (4.1)-(4.2) is an entire func¬ 
tion admitting the following representation 


A(A) = Ao(A)+ /'yi(t)H'i^*dt+ y2(i)C'=^*dt, 

Jo Jo 


(4.8) 


with yi,y 2 G ^^[0,1]. 
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Proof. Consider representations (3.3)-(3.6) for ipjk{-,X), j,k G {1,2}. By Proposition 3.1, 
€ Xi^o(^) IX Xoo,o(Xi), j,k G {1,2}. Therefore by Lemma 2.2, the trace functions 
■) are well defined and summable, -R}'}.(1, •) € T^[0,1], j, k G {1, 2}. Therefore one can sub¬ 
stitute X = 1 in formulas (3.3)"(3.6) and obtain special representations for j,k G {1,2}. 

For instance, 

Pjk{X) := ipjk{'i-,X) = j + j j ^ k. (4.9) 

Inserting these expressions and similar expressions for in (4.5) and taking formula (4.6) 

for Ao(-) into account we arrive at (4.8) with gj{-), j G {1, 2}, being a linear combination of the 
functions i?^(l, ■), j, k G {1, 2}. □ 

In the sequel we need the following definitions (cf. [20]). 

Definition 4.2. (i) A sequence A ;= {Xn}nez of complex numbers is said to be separated if 
for some positive 5 > 0, 

“ -^fcj > 2(5 whenever j / k. (4-10) 

In particular, all entries of a separated sequenee are distinet. 

(a) The sequence A is said to be asymptotically separated if for some no € N the subse¬ 
quence A^q := {Ari}|„|>n(, is separated. 

(Hi) Let A lie in the strip II/i. It is called incondensable if for some L > 0 and N G N 
every reetangle \t — L,t + L] x [—/i, /i] C C eontains at most N entries of the sequenee, i.e. for 
eaeh t gM. the number of integers {n gT, : [Re ^ L, jlm A,i| ^ h} does not exceed N. 

We need the following simple property of incondensable sequences. 

Lemma 4.3. Let A = {Xn}nez be an ineondensable sequence lying in the strip II/j. Then there 
exists eo > 0 and A}) G N such that for any e G (0, eo) every connected component of the union 
of discs U„gzDe(An) has at most Nq discs De(A„). 

Proof. Assume the contrary, i. e., for any e > 0 and K G N there exists connected component 
of He ;= U„gz®e(An) that has at least K discs De(An)- By definition of incondensable sequence 
for some L,N > 0 every rectangle [t — L,t L] x [—h,h], t G M, contains at most N entries 
of the sequence A. Let K > N he some positive integer and pick e to be such that 2Ke < L. 
Consider some connected component of that has M > K discs De(A,i), denote it by C. Let’s 
pick one of the discs Dq = De(A„o), no G Z, in C, and let to = ReA„o. Due to above, rectangle 
[to — L, to + L] X [—ti, h] contains at most N entries of the sequence A. Consider the sequence 
B of all discs Be(A„) in C that have graph distance at most K from Dq. Let’s shows that B 
has more than N discs. If no disc in C has graph distance at least K from Dq, then B contains 
all discs from C and thus cardinality oi B is M > K > N. Otherwise, B has some disc D with 
distance K from Dq. All discs on the path from Dq to D belong to B and hence B has at least 
K > N discs. For each disc ©^(A,!) in B since graph distance from it to Dq is at most K and 
disc radii are e we have |A„ — A„o| < 2Ke. Thus, A„ G [to — L,tQ + L] x [—h, h] since 2K£ < L. 
Thus centers of all discs in B lie in [to — L, to + L] x [—h, h]. Since there more than N discs in 
B it contradicts incondensability property of the sequence A. □ 

To get the asymptotic behavior of the eigenvalues of the problem (4.1)-(4.2) with regular 
boundary conditions we also need the following definition. 
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Definition 4.4. [24] An entire function F(-) of exponential type is said to be of sine-type if 

(i) all zeros of F{-) lie in the strip for some h> t), and 

(ii) there exists Ci,C 2 > 0 and ho > h sueh that 

0 < Cl ^ \F{x + iho)\ ^ C 2 < 00 , X G M. (4-11) 

This definition is borrowed from [24] (see also [20]). It differs from that contained in [25]. 
Namely, it is assumed in [25] that the sequence of zeros of F{-) is separated and the indicator 
function /if(‘) of F{-), 

_ In \F I 

hp{Lp) := lim — - -(/?g(— 7 r, 7 r], (4.12) 

r ^+00 r 

satisfies the condition /ii 7 ’( 7 r/ 2 ) = hp{—'K/2). The latter is imposed for convenience and can 
easily be achieved by multiplication of F[-) by a function with an appropriate 7 G M. 
Recall also the definition of regular boundary conditions. 

Definition 4.5. Boundary conditions (4.2) are called regular if 

JiiJm / 0. (4.13) 

In the case of regular boundary conditions the characteristic determinant Ao(-) has certain 
important properties. 

Proposition 4.6. Let the boundary conditions (4.2) be regular and let A(-) be the characteristic 
determinant of the problem (4.1)-(4.2) given by (4.5). Then the following hold: 

(i) The characteristic determinant A(-) is a sine-type function with /iA(vr/ 2 ) = —bi and 
/ia(— vr/2) = 62 - In particular, A(-) has infinitely many zeros 

A := {Anjnez (4-14) 

counting multiplicities and A C II/i for some /i > 0. 

(ii) The sequence A is incondensable. 

(Hi) For any e > 0 the determinant A(-) admits the following estimate from below 

|A(A)| ^C,(e-'A“Ve-'A“^), AgC\|Jd,(AO, (4.15) 

with some C^ > 0 . 

(iv) The sequence A can be ordered in such a way that the following asymptotical formula 
holds 

, 27m , , ,, - , 

An = 7 - —(1 + 0 ( 1 )) as n —>■ 00 . (4.16) 

02 — bi 

Proof, (i). Let Ao(-) be the characteristic determinant of the problem (4.1)“(4.2) with Q = 0. 
It easily follows from (4.6) that Ao(') admits a representation 

Ao(A) = [ ' F^^daoit), A G C, (4.17) 

Jbi 

with a piecewise constant function o'o(') having precisely four jump-points { 0 , bi,bi + b 2 , 62 }- In 
particular, 


(To (61 + 0 ) - 0 - 0 ( 61 ) = J 32 / 0 and 0 - 0 ( 62 ) - 0 - 0(62 - 0 ) = Ju / 0 . 


(4.18) 
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Let us set 

XrMr,)^ ^£[ 0 , 62 ], 

(4.19) 

and 




a{t) = ao{t) + / g{s)ds. 

(4.20) 


bi 


Combining these notations with formulas (4.8) and (4.17) we arrive at the following representa¬ 
tion for the characteristic determinant 



AgC, 


(4.21) 


It follows from (4.20) and (4.18) that 


cr( 6 i-I-0) - (t(6i) = J 32 / 0 and < 7 ( 62 ) - < 7(62 - 0) = J 14 / 0. (4.22) 


Due to the property (4.22) representation (4.21) ensures that A(-) is a sine-type function with 
/iAo(vr/2) = —bi and h^^{—T:/2) = 62 (see [25]). Moreover, statement (i) is also implied by the 
representation (4.21) (see [23, Chapter 1.4.3]). 

(ii) and (iii). These statements coincide with the corresponding statements of [20, Lemmas 
3 and 4] for sine-type functions (see also [25, Lemma 22.1] in connection with part (iii))- 

(iv) The determinant A(-) belongs to the class A since its zeros lie in the strip II/i (this fact is 
also immediate from representation (4.21)). Therefore it follows from [23, Theorem 1.4.6] that 
for any e G ( 0 , 7 r/ 2 ) 


lim 

t—>-oo 


n±\t) 

t 


27r 

62 - ’ 


(4.23) 


Here n'^\t) = card {n G Z : jAnj < t, I arg An] < e} is the number of zeros of A(-) in 

the domain {z : [arga;] < £,\z\ < t} counting multiplicity, and n^^(t) = 

card{n G Z : |A„| < t, [tt — argA^j < e}. Since A lies in the strip H/i, asymptotic formula (4.16) 

directly follows from (4.23) (see e.g. [48, Proposition 13.1]). □ 


Clearly, the conclusions of Proposition 4.6 are valid for the perturbation determinant Ao(-) 
given by (4.6). Let Aq = {A^j^ez be the sequence of its zeros counting multiplicity. Let us 
order the sequence Aq in a (possibly non-unique) way such that ReA® < ReA®^]^, n G Z. 

Proposition 4.7. Let Q G L^[0,1](8>C^^^, let boundary conditions (4.2) be regular, and let A(-) 
be the corresponding characteristic determinant. Then the sequence A = {Xn}nez of its zeros 
can be ordered in such a way that the following asymptotic formula holds 

An = A°-|- 0 ( 1 ), as n —>■ 00 , n G Z. (4.24) 


Proof. Let s G (0,1). By Proposition 4.6(iii) there exists Ce > 0 such that the estimate (4.15) 
holds. Combining Lemma 4.1 with Lemma 3.5 yields the following estimate 

lA(A)-Ao(A)| < |A| ^ M,. (4.25) 

with certain Mg > 0. Here in the last inequality we have used that bi < 0 < 62 . 

Due to estimates (4.15) and (4.25), the Rouche theorem implies that all zeros of A(-) lie in 
the domain ^ 

Dg := Bm,(0) U Dg, Dg := |J Bg(A°), (4.26) 
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(5.1) 


and in each connected component of the functions A(-) and Ao(-) have the same number of 
zeros counting multiplicity. Since in accordance with Proposition 4.6(ii), the sequence of zeros 
is incondensable, Lemma 4.3 implies that for e small enough each connected component 
of rig contains at most Nq discs D£(A°) with Nq not depending on e. Hence the diameter of 
each connected component of fig does not exceed 2eNQ. Since e > 0 is arbitrary small, the latter 
implies the desired asymptotic formula (4.24). □ 

5. Strictly regular boundary conditions 

Assuming boundary conditions (4.2) to be regular, let us rewrite them in a more convenient 
form. Since J 14 7 ^ 0, the inverse matrix exists. Therefore writing down boundary condi¬ 
tions (4.2) as the vector equation = 0 and multiplying it by the matrix we transform 

them as follows 

iUi{y) = yi{0)+ by 2 { 0 ) + ayi{l) =0, 

\U 2 iy) = dy 2 { 0 ) + cyi{l) + y 2 {l) =0, 

with some a, b,c,d€ C. Now J 14 = 1 and the boundary conditions (5.1) are regular if and only 
if J 32 = ad — be ^ 0. So, the characteristic determinants Ao(-) and A(-) take the form 

Ao(A) = d + + {ad - 6 c)e*^i^ + (5.2) 

A(A) = d + -|- {ad — bc)ipii{X) + ¥? 22 (A) -|- c<y 9 i 2 (A) -|- bip 2 i{X). (5.3) 

Now we are ready to introduce a notion of strictly regular boundary conditions. 

Definition 5.1. Boundary conditions (4.2) are called strictly regular, if they are regular, 
i.e. J 14 J 32 / 0 , and the sequence of zeros Aq = {X^}nez of the characteristic determinant 
Ao(-) is asymptotically separated. In particular, there exists no such that zeros are 

geometrieally and algebraically simple. 

It follows from Proposition 4.7 that the sequence A = {Xn}nez of zeros of A(-) is asymptoti¬ 
cally separated if the boundary conditions are strictly regular. 

Remark 5.2. Let us list some types of strictly regular boundary conditions (5.1). In all of 
these cases the set of zeros of Aq is a union of finite number of arithmetic progressions. 

(i) Separated boundary conditions (a = d = 0, bc^ 0) are always strietly regular. 

(a) Let 61/62 S Q, i-o. 61 = —nib, 62 = n 2 b, ni,n 2 £ N, 6 > 0. Since ad 7 ^ he, Ao(-) is 
a polynomial at of degree ni -\- n 2 . Hence, boundary conditions (5.1) are strictly regular 
if and only if this polynomial does not have multiple roots. In particular, regular boundary 
conditions (5.1) for Dirac operator are strictly regular if and only if {a — df^ 7 ^ —46c. 

Lemma 5.3. Let 6 c = 0 and boundary eonditions (5.1) are regular (i.e. ad 7 ^ 0/. 

(i) Let 

61 In \d\ + 62 In |a| 7 ^ 0. (5-4) 

Then eonditions (5.1) are strictly regular. 

(a) Let 61/62 0 Q. Then eondition (5.4) is necessary for the striet regularity of boundary 
conditions (5.1). 

(iii) Let 61/62 £ Q and condition (5.4) is violated. Then boundary conditions (5.1) are strictly 
regular if and only if 

61 arg(-d) -I- 62 arg(-a) 


271 gcd( 6 i, 62 ) 




(5.5) 
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where gcd( 6 i, 62 ) is the greatest common divisor of real numbers hi, 62 ; ^.e. the largest number 
b > 0 such that bi/b and b^jb are integers. 

Proof, (i) Since be = 0 , the characteristic determinant Ao(-) in (5.2) becomes 

Ao(A) = d + ^ ^ ^ ^ (5.6) 

Let Ai = {Xi^n}nGi. and A 2 = {A 2 ,n}nGZ be the sequences of zeros of the first and second factor, 
respectively. Clearly, 


arg(—a ^) + 27rn .ln|a| 
Ai,n =-;-h i- 


arg(—d) + 27rn .ln|d| 

'^2,n — ^ 


n EZ. 


(5.7) 


bi bi ’ 62 &2 

Thus, Ai and A 2 are algebraically simple and constitute two arithmetic progressions that lie on 
two lines parallel to the real axis. Condition (5.4) written in the form implies that 

these horizontal lines are different. It follows that the sequence of zeros of Ao(-) is separated 
and hence boundary conditions (5.1) are strictly regular. 

(ii) Now assume that a := —bijb^ 0 Q and condition (5.4) is violated. In this case ImAi^„ = 

for n,m € Z, i.e. the progressions Ai and A 2 lie on the same line 


Tm X — — 


In I 


b2 


parallel to the real axis. Hence 

|Ai,n - A2,m| = 2TTbf^\r + n + am\ 


r = 


arg(—a ^) + aarg(—d) 


G 


(5.8) 


Since a is irrational, the Kronecker theorem ensures that for any e > 0 and M > 0 there exist 
n,m such that |n|, \m\ > M and |r + n + am\ < e. This means that the zeros of Ao(-) are 
not asymptotically separated, which proves the result. 

(iii) Finally, assume that condition (5.4) is violated while ol = —61/62 G Q. Now as on the 
previous step, the progressions Ai and A 2 lie on the same line parallel to the real axis and 
condition (5.8) holds. Since 61/62 is rational, the union of the arithmetic progressions Ai and 
A 2 is asymptotically separated if and only if they have no common entries. Due to (5.8) this is 
equivalent to the fact that Diophantine equation n + am = —r does not have integer solutions 
n,m. It is well-known that such equation has solutions if and only if r/gcd(a, 1) G Z. Since 
arg(—a“^) = — arg(—a), this is equivalent to the condition opposite to (5.5), which completes 
the proof. □ 

Remark 5 . 4 . Consider the case 6 = c = 0. It includes periodic (a = d = —1) and antiperiodic 
(a = d = 1) boundary conditions. So, it follows from the statements (i) {see (5.4)) and (ii) that 
the periodic and antiperiodic BC are strictly regular if and only i/ 61 + 62 7 ^ 0. This fact 
demonstrates a substantial difference between Dirac and Dirac type operators. 

The following result demonstrates that in the case 61/62 ^ Q the problem of strict regularity 
of boundary conditions (5.1) is much more complicated than the one discussed in Remark 5.2 
and Lemma 5.3. 

Proposition 5 . 5 . (i) Let a := — 61/62 0 Q, a = 0, 6 c, d G M \ {0}, then boundary condi¬ 
tions (5.1) are strictly regular if and only if 

d 7 ^ —(a-|-1) (|6c|q!““) “+i . (5.9) 

(ii) Let a = 0, be 0, bi = —nff, 62 = 71 . 26 , ni,n 2 G N, 6 > 0 and gcd(ni,n 2 ) = 1. Then 
boundary conditions (5.1) are strictly regular if and only if 

(m + n 2 )”i+" 2 (- 6 c)”L 


(5.10) 
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Proof, (i) It follows from (5.2) that 

Ao(A) = (i-6c-e*^i^ + e*''2A_ (5 1^) 

Let 62 A =; TTx + iy, x,?/ G M. Then e*^ 2 A _ ^nrx-y^ Since bc,d € M, equation 

Ao(A) = 0 is equivalent to the system 


e ^ cos TTx = be ■ e“^ cos avrx — d, 
e~y sin TTX = be ■ e“^ sin avrx. 


From the second equation in (5.12) we have 


In 


y = 


' -bc-smaTTX ^ 
v sin nx } 

a + 1 


be ■ sin avrx 


< 0 . 


Sm TTX 


Substituting it to the first equation we get 


—be ■ sinaTTxA “+i 


sm TTX 

which is equivalent to 


cos TTX = be 


—be ■ sinaTTxA “+i 


sm TTX 


cos avrx — d, 


—6c • sin avrxA sin(Q; + l) 7 rx 


sm TTX 


= -d. 


sm avrx 


(5.12) 


(5.13) 


(5.14) 


(5.15) 


For simplicity we assume that 6 c = —1 and (0 <)a < 1. Then summarizing all previous formulas 
we see that A is zero of Ao(-) if and only if 


62 A = TTX + iy, 

sin aTTX rj 
sinTTir ’ 

ln( 

y = —^ 


x,y e 


sin OLTTX 


a-\-l 


f(x) ■= ('sinovrx'j a+1 sin(Q:+l)7r3; _ 

\ J ’ V sinTTiE / sinoiTra: 


= -d. 


Note that the second relation in (5.16) is equivalent to 


“-((yA-")n(u(“^))) 

u((uo-.....)n(u(5^.|))). 


(5.16) 


(5.17) 


Let us describe the set ^ in a more explicit way. Let n G Z be hxed and let’s find intersection 
A r\ {n,n + 1). Since a < 1, then at most two intervals of the form {mja, {m + l)/a), m G 
Z, intersect with (n, n + 1). There are two cases possible. First, for some m G Z interval 
(m/a, (m + l)/a) fully covers interval (n,n + 1) (i.e. m/a < n < n + 1 < (m + l)/a). In 
this case, if n and m are of the same parity then (n, n + 1 ) ^ A, otherwise ^ H (n, n + 1) = 0. 
The second case, is when for some m G Z we have n<m/a<n + l. In this case it is clear 
that exactly one of the intervals (re, m/a) and {m/a, re + 1) belongs to A depending on parity of 
m — n. Thus, the domain of the 4th equation in (5.16) is the union Dnezin, where /„ is possibly 
empty subinterval of (re, re + 1 ). 
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Put r[x) := sinavrx/sinvrx. Straightforward calculation shows that for x € In-, n € Z, we 
have 


/'(^) = 


ra + l [x] 


{a + 1) sinaTra: 

—t-s+t (x) 

(a + 1) sinavrx 


((a + 1)^ cos(a + l)7rx — (cot ttx + cot aiTx) sin(a + l)7ra:) 

/ \ 2N 


(a + 1)^ sin TTX sin ottx + \/ r(x) cos ttx + 


a cos aixx 




. (5.18) 


Since sinvra: and sinavrx are of the same sign on each interval In it is clear that f'{x) has fixed 
sign on In- Hence, f{x) is strictly monotonic on In- In particular, the equation f{x) = —d, 
X € In, has at most one solution. Denote it by Xn if it exists and let Un be the corresponding 
value of y from the third equation in (5.16). Clearly, all Xn are different and hence all zeros of 
Ao(-) are simple. 

Now let (5.9) be satisfied and assume that boundary conditions are not strictly regular. Since 
zeros of Ao(-) are simple it means that there exists infinite set S' C N such that Xn-i,Xn exist 
for n € S and 

Xn — Xn-i ^ 0, yn — J/n-i —t 0 as n —>■ oo, n€S. (5.19) 

Since zeros of Ao(-) lie in the strip it follows that ly^l < H, n ^ S, with some H > 0. Hence, 
the third relation in (5.16) implies that 

0 < Cl < := < Ca, n G S, (5.20) 

Sm TTXn 

with some C'i,C '2 > 0. Since tn = it is clear that yn — Vn-i ^ 0 is equivalent to 

tn tn—1 ^ O' 

Taking into account the form of the set A described after the formula (5.17) we see that there 
exists unique m = rUn G Z such that either Xn-i < mla < n < Xn oi Xn-i < n < m/a < Xn- 
Moreover, has the same parity as n. Let Si(S 2 ) be the set of those n in S for which the 
first (the second) inequality is satisfied. Since S is infinite and S = Si U S 2 , either Si or S 2 is 
infinite. First consider the case when Si is infinite. For n G Si we put 


don -— mja Xn—l, ^n — ^ mja, din - — Xn Xi- 


(5.21) 


Since Xn-i < mja < n < Xn, then en,<5on,5in > 0. Further, since Xn — Xn-i ^ 0 as n ^ 00 , 
n G Si, then 0, don —> 0, din —?> 0 as n —>■ 00 , n G Si. Hence, for large n G Si we have 

taking into account that m and n are of the same parity 


sinaTTXn sin 7r(m + a(din + Sn)) sm7Ta(din + £n) 
sinvrxn sin7r(n + (5i„) sinvrcJi^ 

sina7ra:„_i sin7r(m — aJon) sinyraJon 

sin TTXn-i sin 7r(n - don - Sn) sin7r((5on + £n} 


(5.22) 

(5.23) 


Here we used the inequality sin art > asinu for 0 < a < 1 and 0 < v < u < tt. 

Since tn — tn-i —> 0, it follows from (5.22) and (5.23) that tn —>■ a, tn-i —>■ a as n —> 00 , 
n G Si. This implies that Snjdin ^ 0 as n —>■ 00 , n G Si. Indeed, since —>■ 0 and din —>■ 0 as 
n —)■ cx), n G Si, then 


a = lim tn = lim 

n—Voo n —^00 

n^S\ n^S\ 


sin7rQ;((5i^ + g^) 
sin irdin 


lim 

n—>-00 
n^S\ 


“1“ ^n) 
din 


a + a lim 

71 —>-00 

n^S\ 


din 


(5.24) 
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Further, note that 

sin(a + l)7rxn sin7r(n + + m + a(5i„ + e„)) sin7r((a + l)(5i„ + 

sm. aiTXn sin7r(m + a(5i„ + £„)) sin7r(Q:5i„ + ac„)) 


(5.25) 


Finally, taking into account the last relation and the fact that £nl^in 0 and ^ a as n ^ oo, 
n G S'!, we have 


—d = lim f(xri) = lim t< 

n—^oo ^ ' n—^oo 

n^Si n^Si 


^ sin(a + l)7rxn ^ i- sin7r((a + l)(5i„ + «£„) 


= a“+i 


sin airxn 

^ l;™ Ol + l+a£n/din _1 « + ! 


lim 

—Voo 
6 Si 


= q;“+i lim 

n—¥oo 

n^Si 


a + aenjdxr, 


= a“+i ■ 


a 


sin7r(Q:5ii 

_ ct 

= {a + l)a “+i. 




(5.26) 


Since be = —1 this contradicts to (5.9). Therefore, zeros of Ao(-) are asymptotically separated. 
The case of infinite S 2 is considered similarly. 

Now let’s prove that opposite statement. As above, for simplicity we assume that be = —1 and 
a < 1. We need to prove that ifd=—(a + l)a “+i then zeros of Ao(-) are not asymptotically 
separated. For n G Z we set m = = [anj. Let S' C Z be some infinite set such that 

rrin — n—)-0asn—>-oo and is of the same parity as n. Since a 0 Q it is clear that such 
set exists. Let us prove that for large enough n G S the equation f{x) = —d has zeros Xn-i-,Xn 
such that Xn-i < m/a < n < Xn and Xn — Xn-i ^ 0 as n ^ 00 . Put 


g{u,v) 


sm au 


1 

“+i sin(au + v) 


smu 


sin au 


(5.27) 


Clearly 


f{Xn-l) = g {TT6on,'^{Son + £n)) , and /(x^) = 5 ( 7 r(( 5 i„ + e^), 7 r( 5 i„) , 

where £n,Son,Sin are defined in (5.21). Due to the special form of d, equation g{u,v) 
equivalent to 


f{u,v) 


sin(ari + v) 
(X \ 


a+l 

— sin V • 



It is easy to prove that 


/(u, 0) > 0,f{u,u) < 0,/(tt,2u) >0, 0 <u < 1/2. 


(5.28) 
= —d is 


(5.29) 

(5.30) 


Hence for each u G (0,1/2) there exists v+ G {u,2u) and V- G (0, u) such that f{u,v±) = 0. 
Clearly, v± = u±(u) is continuous at u. Hence for sufficiently small e > 0 there exists uf,vf 
such that uf — vf = ±e, g{uf,vf) = —d and —>■ 0 as e —>■ 0. Applying this fact for 

shows existence of needed Xn-i and x„. It is easy to prove from f{xn) = f{xn-i) = —d that 
tn ^ a and tn-i — )■ a as n ^ oo where tn-i are defined in (5.22)-(5.23). Which implies 
Un — Vn-i oo and shows that corresponding zeros of Ao(-) are not asymptotically separated. 

(ii) Since 61/62 G Q the set of zeros of Ao(-) is the union of finite number of arithmetic 
progressions. Hence zeros are asymptotically separated if and only if Ao(-) does not have multiple 
zeros, which is equivalent to the fact that Ao(-) and Aq(-) have no common zeros. When a = 0 
wehaveAQ(A) = 7626 *^^''’—f 6 i• 6 c-e*^i'’’. Hence zeros of Ag(-) can be found explicitly. Substituting 
these values into Ao(A) and performing straightforward calculations we see that Ao(-) and Aq(-) 
have no common zeros if and only if condition (5.10) is satisfied. □ 
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6. The Riesz basis property of root vectors system 
6.1. Some auxiliary results. Recall the following definition. 

Definition 6.1. (i) Let 9) be a separable Hilbert space. The vectors system {fk}kez C 9) is 
called Besselian in 9) if {{f, fk)}kez G 

(a) The vectors system {fk} is called a Riesz basis in 9) if it constitutes a basis equivalent 
to an orthonormal one, i.e. there exists a linear homeomorphism T in 9) for which {Tfk} is an 
orthonormal basis. 

(Hi) A sequence of subspaces is called a Riesz basis of subspaces in 9) if there 

exists a complete sequence of mutually orthogonal subspaces and a bounded operator T 

in 9j with bounded inverse such that 9)k = T9)}, A: G N. 

(iv) It is said that a sequence {fk}‘^i of vectors in 9) forms a Riesz basis with parentheses 
if each its finite subsequence is linearly independent, and there exists an increasing sequence 
{nfcl^o ^ ^ such that no = 1 and the sequence 9)k := span{/j}J^~^ ^ constitutes a Riesz basis 
of subspaces in 9). Subspaces 9)k are called blocks. 

Our investigation of the Riesz basis property of the root vectors system of the operator L{Q) 
is heavily relied on the following well-known Bari criterion. 


Theorem 6.2. [17, Theorem VI.2.1] Let 9) be a separable Hilbert space. The vectors system 
{/fcjfcez V 9) forms a Riesz basis in 9) if and only if it is complete, minimal and Besselian in 9), 
and the corresponding biorthogonal system {gk}ke'Z is also complete and Besselian. 

It is well-known that the root vectors system of the operator L{Q) after proper normalization 
is biorthogonal to the root vectors system of the adjoint operator L{Q)*. In this connection we 
give the explicit form of the operator L[Q)* in the case of boundary conditions (5.1). 


Lemma 6.3. Let L{Q) be an operator corresponding to the problem (4.1), (5.1). Then the 
adjoint operator L* := L{Q,Ui,Lf 2 )* is L* = T((5*, R* 2 )); i-a. it is given by the differential 

0 


expression (4.1) with Q*{x) = 


Qi2{x) 


Q2i{x] 

0 


instead of Q and the boundary conditions 


U*i{y) = kbyi{0) + y 2 { 0 ) + dy 2 {l) = 0, 
U* 2 iy) = ayi{0) + yi{l) + k~^cy 2 {l) = 0, 


where k ;= —62^;^^. Moreover, boundary conditions (6.1) are regular (strictly regular) simulta¬ 
neously with boundary conditions (5.1). 

The following lemma plays the key role in the proof of Theorem 1.1. 

Lemma 6.4. Let Q € L^[0,1] 0 and let $(-,A) and 'k(-,A) be the fundamental matrix 
solutions of the equations (4.1) and (3.14) satisfying 4>(0, A) = 'I'(0, A) = I 2 , given by formulas 
(3.2) and (3.15), respectively. Let also $j(-,A) and 'I'j(-,A), j G {1)2}, be the columns of 
these matrices {cf. (3.2) and (3.15)). Then for any incondensible sequence {pn}nez the systems 
{4>j(-,/r„)}„6Z and are Besselian in ^^[0, 1] (gjC^, j G {1,2}. 
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Proof. Consider the case of the system /rn)}nGZ- Let / := col(/i,/ 2 ) S L^[0,1] (8) C^. 

Taking into account formulas (3.3), (3.5) we get 

(/,^>l(-,/rn))L2[0,l]<^C2 = [flix)(pu{x, Hn) + f2{x)<f2lix, dx 

= J /i(a;)e*^i/^'‘*(ix + ^ J fj{x)^J dx, ( 6 . 2 ) 

where := i2+i(-,-), Nj 2 {-,-) ■= Rj 2 {-r)i 3 ^ {1^2}. Further, note 

J fjix)^J A(jfc(x, dx = J gjk{t)e'^^kr-nt ( 6 . 3 ^ 


where 


9jk{t) ■= Njk{x,t)fj{x)dx, j,ke {1,2}. 


(6.4) 


By Proposition 3.1, Njk{-,-) S Xi(n) nXoo(n). Therefore by Lemma 2.1 the Volterra type 
operators 

Njk-f^ [ Njk{x,t)f(t)dt and N*f. : f ^ f Njk{x,t)fj{x)dx 
Jo Jt 

are bounded in L^[0,1], hence gjk G L^[0,1], j,k G {1,2}. Taking this inclusion into account 
and inserting expressions (6.3) into (6.2) one rewrites this equality as 


(/,‘^’l(‘,/^n))2,2[ou](g)C2 — e 


iblflnX 


L2[0,1] 




j,k=l 




L2[0,1] 


(6.5) 


Since the sequence {/injnez is incondensible, then by [20, Lemma 2 ] the sequence of exponents 
is Besselian in L^[0,1]. The latter implies 

{(ffife,e*^'“^"*)L2[o^i]}nez G /^(Z) and {(/,e*^''^"*)L2[o,i]}nez G /^(Z). 

Combining these inclusions with representation (6.5) shows that the system {<l)i(-,/r„)}„gz is 
Besselian. The systems of functions {< 1 > 2 (-,//n)}nez and j{-,jm)}n&'L are treated similarly. □ 


6.2. Proof of the main result. Now we are ready to prove the main result of the paper. 

Proof of Theorem 1.1. According to Proposition 4.7 and Definition 5.1, the operator L{Q) has 
countably many eigenvalues {\n}nez- Moreover, they are of finite multiplicity, asymptotically 
simple and separated, and are located in the strip Ll/i = {A G C : jlm A| ^ h}. 

Since boundary conditions are regular, one can transform them to the form (5.1). 

(i) In this step assume that | 6 | + jcj 7 ^ 0. Without loss of generality it suffices to consider the 
case 6/0. Let ^ = {fn}nez and <5 = {gn}nez be the system of root vectors of the operators 
L{Q) = L{Q,Ui,U 2 ) and L{Q)* = L(Q*, ?7*i, 17*2)), respectively. By [34, Theorem 1.2], each of 
the systems ^ and © is complete and minimal in L^[0,1] ( 8 > C^. Therefore these systems can be 
chosen to be biorthogonal to each other. 

First we indicate the explicit form of the functions fn and gn for n large enough. 
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To this end one easily gets from (3.2) and (5.1) that 

f/i(^>i(-, A)) = (/?ii(0, A) + bip2i{0, A) + a(/?ii(l, A) = 1 + a(/?ii(A), ( 6 . 6 ) 

Ui{^2{-, A)) = (/?i2(0, A) + bip22{0, A) + aifuil, X) =b + aipuiX). (6.7) 

Since ifuiXn) = o(l) as n —>■ oo (see (3.12)) and 6/0, one gets that b+aifuiXn) / 0 for |n| > ni 
with some rei € N. Therefore one derives from ( 6 . 6 ), (6.7), and (3.2) that the vector-function 


fn {-) :=u, ($2(-, A„)) $i(-, A„) - Ui ($i(-, A„)) $2(-, A„) 

=(6-Fa(/3i2(An))^>i(-,An) - (l + a(/?ii(A„))^> 2 (-,A„), |n| > m, (6.8) 

is a non-trivial eigenfunction of the operator L{Q) corresponding to the eigenvalue A^. Since 
boundary conditions (5.1) are strictly regular, it follows from Proposition 4.7 that the sequence 
A = {Xn}n£Z of the eigenvalues of L{Q), i.e. the zeros of A(-), is asymptotically separated. 
In particular, there exists no G N such that the eigenvalues of L{Q) are geometrically and 
algebraically simple. Therefore fn is the unique up to a multiplicative constant eigenfunction of 
the operator L{Q) corresponding to A„ for |n| > max{no,ni}. 

Similarly, one easily gets from (6.1) that 

C/,i(Ti(-,A)) = A:6 + #2 i(A), C/*i(T2(-, A)) = 1 + #22(A). (6.9) 

Moreover, in accordance with (3.12) there exists n 2 G N such that kb+d'ip 2 i{Xn) / 0 for |n| > n 2 . 
Therefore Lemma 6.3 ensures that the vector-function 


5n(-) := (^2(-, A„))Ti(-, A„) - [7*1 (T/-, A„))T2(-, A„) 

= ( 1 d'i/’22(An))^i(-, A„) - (fc6-hd/2i(An))^2(-, A„), |n| > n2, ( 6 . 10 ) 

is a non-trivial eigenfunction of the operator L[Q)* corresponding to the eigenvalue A„. Since 
the eigenvalues of L{Q)* constitutes a sequence {Xn}nez, they are geometrically and alge¬ 
braically simple simultaneously with {Xn}nez, he. for |n| > no. Therefore gn is the unique 
up to a multiplicative constant eigenfunction of the operator L{Q)* corresponding to An for 
|n| > max{no, n 2 }. 

Further, it follows from (3.12) and (3.17) that 

ll/nf = |6 + ao{l)f • ||$i(-, An)lP + |1 + + o(l)p • ||4>2(-, An)f 

-2Re (6 + ao(l))(l + + o(l))(4«i(-, An), An)) 

= |6p • ||<I>i(-, An)|P + |1 + ae*^^''‘’*|^||<I> 2 (-, An)|p + o(l) as n —>■ oo. (6-11) 

Similar relation is valid for gn- Combining these relations with estimates (3.18) and noting that 
An G n,,, n G N, yields 

||/n|| X 1, IbnII 1, for large n. (6.12) 

Here the symbol On x 1 for large n means that there exists no G N and Ci,C 2 > 0 such that 
Cl < I On I < C 2 , |n| > no. In particular, vector-functions fn, gn are non-zero for large n. 
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Similarly starting with (6.8), (6.10), and using relations (3.12), (3.16), and noting that fe € R, 
one arrives at the following important asymptotic relation 

{fn,9n) = (6 + a-o(l))(l + de“*^2^" + o(l)) ($i, ^'i) 

+ (l + + o(l)) (kb + d- o(l)) ($2, ^-2) 

= +/c + 1) + o(l) as n ^ oo. (6.13) 


Further, using the formula (4.6) for the perturbation determinant Ao(-) one easily gets 
^ (Ao(A)e-*^i^) = ^ + ad-bc + 

= + k + l). (6.14) 


Since boundary conditions (4.2) and (5.1) are equivalent, the latter are also strictly regular. 
Therefore in accordance with Dehnition 5.1 the sequence of zeros Aq := {A^jnez of the deter¬ 
minant Ao(') is asymptotically separated. In other words, there exist uq a N and d > 0 such 
that the separation condition (4.10) is satished for |j|, |A:| > no- 

Combining estimate (4.15) with the Minimum and Maximum Principles for analytic functions 
(cf. [25, Lemma 22.2]), we arrive at the following two-sided estimate 

Cl < |A[,(A)| <C 2 , AgD 5 (A°), |n|>no. (6.15) 

with certain Ci = Ci(d) > 0, C 2 = 02 ( 6 ) > 0. It follows from asymptotic formula (4.24) that 
An = A° -|- 0 ( 1 ) G ©^(An) for n large enough. Combining this inclusion with estimates (6.15) 
yields AQ(An) x 1 for large n. Since jlm An] ^ h, relations (6.13), (6.14), and (6.15) imply 


(/njffn)-! for large n. 

Thus we can normalize the systems {/n} and {gn} by putting 


fn ■ = 


fn 


ll/nir ' (/n,5n)’ 

Clearly, ||/n|| = 1 and (/n,5n) = 1 for jnj > m, i.e. the sequences 

■ {fn}\n\>m and Q . {9n}\n\>m 

are biorthogonal. It follows from (6.8), (6.10) and (6.17) that 

^ 6 + ac/Pl2(An) l + aipii{\n) , 

JnV) - -||-^-|j-4>l(-,An)--4'2(-,An), 

- ^ ||/n||(l +dV;22(An)) .j,^^ _ ||/n || (^6#21 (An)) 


ll/nllffn 


n| > m := max{no, ni, 712 }. 


(/nj 9 n) 


{fn,9n) 




(6.16) 

(6.17) 

(6.18) 

(6.19) 

( 6 . 20 ) 


By Proposition 4.6, the sequence of eigenvalues A = {Anjnez is condensible. Therefore, by 
Lemma 6.4, the sequences of vector-functions {#(•, An)}|n|>m and {#(•, An)}|n|>m! j ^ {1)2}, 
are Besselian in #[0,1] 0 C^. It follows from (3.12), (6.12), and (6.16) with account of the 
inclusion A C B/i, that the coefficients at #■(•, A^) and 'I'j(-, A^), j G {1,2}, in (6.19) and (6.20) 
are bounded in n. Hence, the sequences J- and Q given by (6.18) are also Besselian in L^[0, 1](8>C^. 

Clearly, the systems 5 and © contain and Q, respectively. Since 5 differs from by a finite 
number (in fact, 2m + 1) entries and J- is Besselian in L^[0,1] 0 C^, the system ^ is Besselian 
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too. Similarly <5 is Besselian in L^[0,1] 0 C^. The Riesz basis property of the system 5 in 
L^[0,1] (8> is now implied by the Bari criterion (Theorem 6.2). 

(ii) Now assume that 6 = c = 0. In this case ad ^ 0 and Ao(A) = {d + 

Let Aj* = {Aj „}nez nnd = {A 2 „}nGZ be the sequences of zeros of the hrst and second 
factor, respectively. Clearly, these sequences constitute arithmetic progressions lying on the lines 
parallel to the real axis. Since the boundary conditions are strictly regular, by the assumption, 
the arithmetic progressions A)* and A 2 are separated, i.e., |A5„ “ A^^I > 25, m, n G Z for some 
5 > 0. This implies the following asymptotic relations 

1 + X 1, d + ^ x for large \n\. (6-21) 

In accordance with Proposition 4.7 the sequence A = {AnInez of zeros of A(-) admits a de¬ 


composition A = Ai U A 2 where Kj = {Xj^n}ne'E, j £ {15 2}, and Xj^n = A°„ + o(l) as n —>■ 00 . 
Combining this representation with Proposition 3.4 and (6.21) yields 

1 + a(pii{Xi^n) - 1, d + (P 22 i>^ 2 ,n) - 1 for large |n|. (6.22) 

Similarly to (6.22) one derives 

o +'!/'ii(Ai,n) - 1, (1 + d^/' 22 (A 2 ,n)) - 1 for large |n|. (6.23) 

Combining this relation with (6.22) yields the existence of n} G N such that 

(1a(^ii(Ai,n))(a + V'ii(Ai,n)) 7 ^ 0 for |n| > n{. (6.24) 

Therefore it follows from (6.8) (with b = 0 and a 0) that for |n| > n{ the vector-function 

/l,n(-) := 0 ¥^ 12 (Al,n)‘hl(-, Ai^n) — (1 -|- (Ai,n))‘h2(', Ai^n), (6.25) 

is the non-trivial eigenfunction of the operator L{Q) corresponding to the eigenvalue Ai^n- More¬ 
over, it follows from the second equation in (6.1) and (6.24) that the vector-function 

9l,n{-) '■= tl*2(^2(-, Ai^„))'I'i(-, Ai^n) — C/* 2 ('I'i(-, Ai_„))'I' 2 (-, Ai^n) 

:= V'i 2 (Ai,n)^i(-, Ai,n) - (a + ^ii(Ai,n))T 2 (-, Ai,„), |n| > n'l, (6.26) 

is a non-trivial eigenfunction of the operator L{Q)* corresponding to the eigenvalue Ai^^. Since 


boundary conditions (5.1) are strictly regular, Proposition 4.7 ensures existence of Uq G N such 
that the eigenvalues {Ai^,„} of L{Q) are geometrically and algebraically simple for |n| > Uq. 
Therefore /i,n(‘) (9i,n(')) is the unique up to a multiplicative constant eigenfunction of the 
operator L{Q) {L{Q)* ) corresponding to the eigenvalue Ai^„ (Ai,n) for \n\ > max{nQ,n{}. 

Further, combining (6.22) with relation (6.23) for large enough \n\ and applying. Proposi¬ 
tion 3.4 and Corollary 3.6, we arrive at the following asymptotic relations 

Il/i,n|| - 1, ||5i,n|| - 1, (/i,n,5i,n) - 1 for large n, (6.27) 

(cf. (6.12) and (6.13)). Performing normalization of the sequences {/2,n}nGZ and {g 2 ,n}ne‘Z in 
the same way as for the sequence (6.17) and repeating the same argument we get that their 
normalizations are besselian biorthogonal sequences. 

Going over to the second branch {A 2 ,n}nGZ of eigenvalues we obtain from (3.2) and (5.1) with 
account of the assumption c = 0 that 

/2,n(-) := U 2 (‘h 2 (-, A2,n)) ‘hi(-, A2,n) “ U 2 (‘hi(-, A2,n)) ‘h 2 (-, A2,n) 

= {d + ‘/^22(A2,n))‘hl(-, A2,n) “ 7^21 (A2,n)‘h2(-, A2,n) 


(6.28) 
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It follows from the second relation in (6.22) that for n big enough f 2 ,n{') is a non-trivial eigen¬ 
function of the operator L{Q) corresponding to the eigenvalue A 2 ,n- Similarly, it follows from 
the first equation in (6.1) (with 6 = 0) and (6.23) that for n big enough the vector-function 

92,n{-) ■= U^i ('I'2(-, A2,n))'I'l(-, )^2,n) “ ('I'l(', A2,n))^2(-j A2,n) 

= (1 -|- dV' 22 (A 2 ,n))'I'l(-j A 2 ,n) “ (A 2 ,n)'I' 2 ('j A 2 ,n) (6.29) 

is a non-trivial eigenfunction of the operator L{Q)* corresponding to the eigenvalue A 2 ,n- Starting 
with (6.28) and (6.29) one completes the proof in the case of strictly regular BC by repeating 
the above reasonings. □ 

6.3. The case of general potential matrix. In applications systems (4.1)-(4.2) appear with 
potential matrices having non-trivial diagonal, i.e. of the form 

Q = G L^[0, 1] 0 C2x2. (6.30) 

First we apply gauge transformation to reduce system (4.1)-(4.2), with a potential matrix Q{-) 
of the form (6.30) to similar system with off-diagonal potential matrix Q. To this end we put 

Wj{x) := exp (^—ibj J Qjj{t)dt^ , x G [ 0 , 1 ], j G { 1 , 2 }, (6.31) 

Lemma 6.5. Let Q he a summable matrix given by (6.30). Then the operator L{Q) = 
L{Q,Ui,U 2 ) is similar to the operator L{Q) = L{Q,Ui,U 2 ) given by (4.1)-(4.2) with the same 
B, a potential off-diagonal matrix Qf), 

, fc(i);=u.r‘W»2{x), (6.32) 

instead of Q and the boundary conditions 

Uj{y) := ajiyi(O) -h aj2y2(0) -h u;i(l)aj 3 yi(l) -h rc2(l)ai4y2(l) = 0, j G {1, 2}. (6.33) 

Proof. See the first part of the proof of [30, Proposition 3.4]. □ 

Corollary 6.6. Let Q he a summable potential matrix given by (6.30) and let boundary condi¬ 
tions (4.2) he separated and regular, i.e. 

ayi(O) -I- 61 / 2 ( 0 ) = cyi{l) -|- dy 2 {l) = 0 and abed / 0. 

Then the eigenvalues of the corresponding operator Lc,d{Q) o,re asymptotically separated and 
there exists /3 G C such that the following asymptotic formula holds 

^ 27rn ^ . . 

A„ = 7 -—-l-p-l-o(l) as jnj 00 . 

62 - 61 

Moreover, the system of root vectors of Lc,d{Q) forms a Riesz basis in L^[ 0 , 1 ] 0 C^. 

Lemma 6.7. Let boundary conditions (4.2) be regular. Then there exists w 0 such that the 
boundary conditions 

Uj{y) := ajiyi{0) + 0 ^ 22 / 2 ( 0 ) + waj 3 yi{l) + 0 ^ 42 / 2 ( 1 ) =0, j € {1, 2}, (6.34) 

are strictly regular. 
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Proof. Since boundary conditions (4.2) are regular, J 14 J 32 7 ^ 0, we can assume without loos 
of generality that J 14 = 1 . Let Ao(-) be the characteristic determinant corresponding to the 
BVP (4.1), (6.34). It is easily seen that 

Ao(A) = J12 + ( 6 . 35 ) 

Let us assume for simplicity that 61/62 € Q, i.e. 61 = —nib, 62 = n 2 b, ni,n 2 G N, 6 > 0. In this 
case we can rewrite Ao(-) in the following form 

Ao(A) = P^iz) := + wJs^z^^ + wJ32- (6.36) 

Since J 32 7 ^ 0, one can choose w 0 such that the zeros of the polynomial Pw{z) are sim¬ 
ple. Clearly, for such w the zeros of Ao(-) are asymptotically separated and thus boundary 
conditions (6.34) are strictly regular. □ 

To state the next result we recall that ma{\o) and mg(Ao) denote the algebraic and geometric 
multiplicities of Aq, respectively. Moreover, if Aq is an isolated eigenvalue, then ma(Ao) equals 
the dimension of the Riesz projection. 

We need the following known abstract result (see e.g. [44]). 

Proposition 6.8. Let L he an operator with compact resolvent in a separable Hilbert space io and 
let {Xn}nez be the sequence of its distinct eigenvalues. Assume that ma{Xn) < 00 for n G N and 
that A has finitely many associative vectors, i.e. there exists no G N such that ma{Xn) = mg{Xn) 
for jnj > no. Further, assume that 

|An| ^ Cjn], [ImAnl^r, n G Z, (6.37) 

for some C, r > 0. Finally, let the system of root vectors of the operator L forms a Riesz basis 
in Sj. Then for any bounded operator T in Sj the system of root vectors of the perturbed operator 
A = L + T forms a Riesz basis with parentheses in i^. 

Proof. Since L has hnitely many associated vectors, there exists a finite-dimensional operator 
K such that the operator L + K has no associative vectors, i.e. nia{Xn) = rng{Xn) for n G N. 
Then the system of eigenvectors of L-|-iL constitutes a Riesz basis in f), i.e. the operator L + K 
is similar to a normal operator H. The latter admits a representation H = Hn iHj where the 
operators Hn := [H + H*)/2 = and Hi := {H — H*)/2i are self-adjoint and commute (see [7, 
Theorem 6.6.1]). Since the spectrum of H lie in a strip 11,-, its imaginary part Hj is bounded, 
||L^/|| A T. Clearly, inequality (6.38) remains valid for eigenvalues of H maybe with another 
constant C > 0. Therefore the operator A is similar to a bounded perturbation of the self-adjoint 
operator H whose eigenvalues satisfy (6.38). Hence, by [19, Theorem 3.1] (Katsnel’son-Markus- 
Matsaev theorem, see also [36, 37]) the system of root vectors of A constitutes a Riesz basis 
with parentheses. □ 

Proposition 6.9. Let Q be a summable potential matrix given by (6.30) and let Lc,d{Q) be 
the operator associated in L^[0,1] ( 8 > with the BVP (4.I)-(4.2). Assume that boundary con¬ 
ditions (4.2) are regular. Then root vectors system of the operator Lc,d{Q) forms a Riesz basis 
with parentheses in Lfif), I] (g) C^. 

Proof. It is clear that the regularity of boundary conditions is preserved under gauge transfor¬ 
mation used in Lemma 6.5. Therefore one can assume that Q is off-diagonal. Now let us consider 
a perturbation of the operator L{Q) by a constant diagonal potential matrix Qo = diag(( 7 o, 0 ), 
go G C. Applying Lemma 6.5 again we see that the operator Lc,d{Q + Qo) is similar to the 
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operator j^{Q) with off-diagonal Q and with boundary conditions C?/(0) -|- Dy{l) = 0, where 
D = D ■ diag(r(;, 1), w = . By Lemma 6.7, we can choose qq € C such that the boundary 

conditions Cy{0) -|- Dy{l) = 0 are strictly regular. Let us verify that the operator L{Q + Qq) 
satisfies conditions of Proposition 6.8. Since eigenvalues of L{Q + Qo) are asymptotically sep¬ 
arated, it follows that their algebraic multiplicities are finite and L(Q -j- Qq) has finitely many 
associated vectors. 

According to Theorem 1.1 the root vectors system of the operator L(Q + Qo) forms a Riesz 
basis in L^[0,1] (8>C^. On the other hand, inequalities (6.38) are implied by Proposition 4.6(iv). 
Thus, the operator L(Q + Qo) satisfies the conditions of Proposition 6.8, and hence the root 
vectors system of the original operator L(Q) = L{Q + Qo) — Qo forms a Riesz basis with 
parentheses. □ 

Remark 6.10. (i) Note that inequalities (6.38) are valid for the roots of any entire function of 
exponential type a < oo with infinitely many zeros. Namely, the following inequalities hold 

Tl 

|An| ^ —, \n\ > N, (6.38) 

ea 

for all but finitely many numbers. In particular, they are valid for the roots of A{-). However, 
Proposition 4.6 (iv) gives sharp asymptotic. 

(a) The Riesz basis property for 2x2 operators Lc,d{Q) with separated boundary conditions 
was established earlier than for the operators with general regular boundary conditions. Namely, 
this property was proved firstly in [55] and later on in [8] and [3] for B = diag(—1,1), Q € 
L^[0,1] ® , and in [18] for B = diag(6i, 62), Q & C'^[0,1] O 

(Hi) The Bari-Markus property of the Riesz projectors of unperturbed and perturbed BVPs 
for separated, periodic and antiperiodic boundary conditions was established in [8] and reproved 
by another method in [3]. In [9] similar results have been obtained for general regular boundary 
conditions. Finally, in the recent paper [40] the results of [8] regarding the Bari-Markus property 
in L^[0,1] O were extended to the case of the Dirichlet BVP for 2m x 2m Dirac equation 
with Q G L2([0,l];C2"^^2m^_ 

7. Application to the Timoshenko beam model 
Consider the following linear system of two coupled hyperbolic equations for t ^ 0 

Ip{x)^tt= Kix){W^-<^) + {EI{x)^^)^-pi{x)<^t, xG[0,^], (7.1) 

pix)Wu= {Kix)iW,,-<!>)),,-p 2 {x)Wt, xe[0,i]. (7.2) 

The vibration of the Timoshenko beam of the length i clamped at the left end is governed by 
the system (7.1)-(7.2) subject to the following boundary conditions for t ^ 0 [53]: 


lT(0,t) = $(0,t) = 

0, 

(7.3) 

{EI{x)^x{x,t) -\-ai^t{x,t) + l3iWt{x,t))\^^i = 

0, 

(7.4) 

{K{x){Wx{x,t) - ^(x,t)) +a2Wt(x,t) -h /92^t(x,t))l^^^ = 

0. 

(7.5) 


Here W{x,t) is the lateral displacement at a point x and time t, ^{x,t) is the bending angle 
at a point x and time t, p{x) is a mass density, K{x) is the shear stiffness of a uniform cross- 
section, Ip{x) is the rotary inertia, EI{x) is the flexural rigidity at a point x, pi{x) and P 2 {x) 
are locally distributed feedback functions, cej,fij G C, j G {1,2}. Boundary conditions at the 
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right end contain as partial cases most of the known boundary conditions if ai,a 2 are allowed 
to be infinity. 

Regarding the coefficients we assume that they satisfy the following general conditions: 

p,Ip,K,EIeClO,£j, pi,p 2 eL^lO,£], (7.6) 

0 < Cl ^ p(x),I^(x),K(x),EI(x) 4 C 2 , xG[0,£]. (7.7) 

The energy space associated with the problem (7.1)~(7.5) is 

f) := H^[0,£] X L'^iOJ] X H^[0,£] x (7.8) 


where Hq[0,£] := {/ G : /(O) = 0}. The norm in the energy space is defined as follows: 


III 


rl 

/ {EI\y'i\^ + Ip\y 2 ? + K\y'.^-yi\^ + p\y 4 \^)dx, y = col(yi,y2,2/3,2/4)- 
7o 


The problem (7.1)-(7.5) can be rewritten as 

yt = iCy, y{x,t)\t=o =yo{x), 

where y and C are given by 


y = 


on the domain 

dom(£) = 


/ $(x,t) \ 


^2/1^ 


^tix,t) 

c 

2/2 

1 

W{x,t) 

5 

2/3 

i 

\Wtix,t)) 


\y^/ 



/ 


V 


2/2 


^(^K{x){y'^-yi) + {EI{x)y[y -pii. 

2/4 

^((^(a:)(2/3 - 2/i))'-P2(a:)2/4) 


/ 


[y = col(yi,y2,2/3,2/4) : 2/i,2/2,2/3,2/4 S , 

EI-y[€AC[0,£], {EI-y[y -piy2€L^[0,£], 
K-{y'^- yi) G AC[£),£], {K ■ (y' - yi))' -p22/4 G L‘^[0,£], 
{El ■y'i){£) + aiy2{£) + / 3 iy 4 (^) = 0 , 
{E ■ (ys - yi)){£) + a2yi{£) + P2y2{£) = o|. 


(7.9) 

(7.10) 


(7.11) 


(7.12) 


Timoshenko beam model is investigated in numerous papers (see [53, 22, 49, 58, 57, 56] and 
the references therein). A number of stability, controllability, and optimization problems were 
studied. Note also that the general model (7.1)-(7.5) of spatially non-homogenous Timoshenko 
beam with both boundary and locally distributed damping covers the cases studied by many 
authors. Geometric properties of the system of root functions of the operator £ play important 
role in investigation of different properties of the problem (7.1)-(7.5). 

Below we establish the Riesz basis property with parentheses of the operator £, without 
analyzing its spectrum. For convenience we impose the following additional algebraic assumption 
on £: 

v{x) := = const, x G [0,£], (7.13) 

Clearly, (7.13) is satisfied whenever Ip{x) = Rp{x), where R = const is a cross-sectional area 
of the beam, El and K are constant functions, while p G AC[£),£] and is arbitrary positive (cf. 
condition (7.19)). Our approach to the spectral properties of the operator L is based on the 
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similarity reduction of £ to a special 4x4 Dirac-type operator. To state the result we need 
some additional preparations. 

Let 7 (-) be given by 


I j (x) 

\ ^ = 6 i 7 (x), where 6 i > 0 and / ^{x)dx = 1 . 

V EI{x) Jo 


(7.14) 


Conditions (7.6) and (7.7) imply together that 7 G C[0,.^] and is positive. Further, in view 
of (7.13) we have 


' p{x) 
K{x) 


= b 2 'y{x), where 62 > 0 . 

diag(- 6 i, 61 ,- 62 , & 2 ). 
-2i diag{Ip{x),Ip{x),p{x),p{x)), 

^EI{x)Ip{x), h 2 {x) := s/Klxjf^. 


Let 

B := 
e{x) := 

hi{x) := 

In the sequel we assume that 

/ii,/i 2 G AC[d 4 ]. 

Therefore, according to (7.6)~(7.7) the following matrix function is well-defined: 

(pi + /i( Pi - h[ h2 -6-2 \ 


Q{x) := 0 (x) 


Pi + h[ Pi - h[ /12 - 6-2 

-/l 2 -/l 2 P2 + 6-2 P2 - h '2 

\ 6,2 /l 2 P2 + h '2 P2 - h' 2 ) 


Next, we set 


nx 

t{x) = / ^{s)ds, xG[0,£]. 

Jo 


(7.15) 

(7.16) 

(7.17) 

(7.18) 

(7.19) 

(7.20) 

(7.21) 


Since 7 G C[d,(\ and is positive, the function t{-) strictly increases on [0,^], t(-) G C^[0,£], 
and due to (7.14) = 1. Hence, the inverse function x(-) := is well defined, strictly 

increasing on [0,1], and x(-) G £*^[0,1]. Next, we put 


Q{t) := Q{x{t)) =: {qjk{t))lk=i, t G [0,1]. 


Finally, let 


C = 


/I 

1 

0 

0\ 


/ 

0 

0 

0 

0 \ 

0 

0 

0 

0 

, ^ = 

ai 

-hi(£) 

ai -L hi{l) 

Pi 

Pi 

0 

0 

1 

1 


0 

0 

0 

0 

VO 

0 

0 

0^ 


\ 

/I2 

h 

Oi2 - 6-2(1') 

02 + 62(1) y 


(7.22) 


(7.23) 


Proposition 7.1. [29, 30, Proposition 6.1] Let functions p, Ip, K, El,pi,p 2 ,hi,h 2 satisfy con¬ 
ditions (7.6), (7.7), (7.13) and (7.19). Then the operator C is similar to the 4x4 Dirac-type 
operator L := Lc,d{Q) with the matrices B,C,D, and Q{-) given by (7.16), (7.23) and (7.22). 


Theorem 7.2. Let conditions (7.6), (7.7), (7.13), (7.19) be satisfied and let also 

/3i =/32 = 0, aif^±hi{£), 02 / ±^2(-^)- (7.24) 

Then the system of root functions of the operator L forms a Riesz basis with parentheses in Sj. 
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Proof. Consider the operator Lc,d{Q) defined in Proposition 7.1. Since /3i = /32 = 0 we can 
represent it as bounded perturbation of the direct sum of two 2x2 Dirac operators: 

LcMQ) = L{UuVi,Qi) (B L{U2,V2,Q2) + Q, (7.25) 

HQj) = -i(^Q^ ^^y' + Qiy, y = coi(yi,?/2), (7.26) 

Ujiy) := yi(0) +y2(0), Vj{y) := (aj - hj{i))yi{l) + {oj + hj{£))y 2 {l), j G {1,2}, 

(7.27) 

Qi{t) = Q,{x{t)), 4 = (-2f/p)-'(^| + J^/ pl-hj)’ 

«.(*)=<3.w*)). &=(-2v)-‘(}j};:? p.29) 

Q(t) = Q(x{t)), § = e-lcodiag((^J;2 (7.30) 

It follows from (7.6), (7.7), (7.19) that <5i,Q2 G -^^[0,1] (g) and Q G L°°[0,1] x 
Due to conditions (7.24), the operator L{Uj,Vj,Qj) is a 2 x 2 Dirac operator with separated 
regular boundary conditions. By Corollary 6.6, the system of its root vectors forms a Riesz basis 
in L^[0,1] (g) and its eigenvalues have a proper asymptotic, in particular, inequality (6.38) 
is satisfied for them. It is also clear that L{Uj,Vj,Qj) has finitely many associated vectors. 
Clearly, the direct sum L := L{Ui, Vi, Qi) 0 L{U 2 , V2, Q 2 ) has the same properties. Since Q is 
bounded, the operator Lc,d{Q) is a bounded perturbation of ’’spectral” operator L. Hence by 
Proposition 6.8, the system of root vectors of the operator Lc,d{Q) forms a Riesz basis with 
parentheses in L^[0,1] 0 C^. Since, by Proposition 7.1, L is similar to the operator Lc^d{Q)-, 
the system of root functions of L forms a Riesz basis with parentheses in i^. □ 

Remark 7.3. In [29, 30] the same result was proved under additional smoothness assumptions 

Pi,P2€L°°[0,i], hi,h2 GLipJO,^]. (7.31) 

Hence Theorem 7.2 is considerable generalization to the most general conditions on coefficients. 
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